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^ ■ Abstract 



X 



Ground-state properties of exotic even-even nuclei with extreme neutron-to- 
proton ratios are described in the framework of the self-consistent mean-field 
5^ I theory with pairing formulated in coordinate space. This theory properly 

accounts for the influence of the particle continuum, which is particularly im- 
portant for weakly bound systems. The pairing properties of nuclei far from 
stability are studied with several interactions emphasizing different aspects, 
such as the range and density dependence of the effective interaction. Mea- 
surable consequences of spatially extended pairing fields are presented, and 
the sensitivity of the theoretical predictions to model details is discussed. 
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I. INTRODUCTION 



One of the most exciting challenges in today's nuclear structure is the physics of exotic 
nuclei far from the line of /3-stability. What makes this subject particularly interesting 
(and difficult) is the unique combination of weak binding and the proximity of the particle 
continuum, both implying the large diffuseness of the nuclear surface and extreme spatial 
dimensions characterizing the outermost nucleons HTj-Q. 

For the weakly bound nuclei the decay channels have to be considered explicitly. Due 
to the virtual scattering of nucleons from bound orbitals to unbound scattering states, the 
traditional shell-model technology becomes inappropriate. The proper tool is the continuum 
shell model f^,^] which correctly accounts for the coupling to resonances; the single-particle 
basis of the continuum shell model consists of both bound and unbound states. The ex- 
plicit coupling between bound states and the continuum and the presence of low-lying \ow-£ 
scattering states invites strong interplay between various aspects of nuclear structure and 
reaction theory. 

Particularly exciting are new phenomena on the neutron-rich side. Because neutrons 
do not carry an electric charge, the neutron drip line is located very far from the valley of 
/3-stability. Consequently, the neutron drip-line systems (i.e., those close to the neutron drip 
line) are characterized by unusually large N/Z ratios. The outer zone of these nuclei are 
expected to constitute essentially a new form of a many-body system: low-density neutron 
matter (neutron halos and skins). 

Except for the lightest nuclei, the bounds of neutron stability are not known exper- 
imentally. Theoretically, because of their sensitivity to various theoretical details (e.g., 
approximations used, parameter values, interactions) predicted drip lines are strongly 
model-dependent. The placement of the one- neutron drip line, defined by the condition 
Sn{Z,N) = Bn{Z , N) — Bn{Z , N — 1) = 0, is solely determined by the binding energy differ- 
ence between two neighboring isotopes. Analogously, the vanishing two-neutron separation 
energy, S2n{Z, N) = Bn{Z, N) — Bn{Z, N — 2), defines the position of the two-neutron drip 
line. Since experimental masses (binding energies) near the neutron drip lines are unknown, 
in order to extrapolate far from stability, the large-scale mass calculations are usually used 
(see, e.g., |17|-[TI|). However, since their techniques and parameters are optimized to reproduce 
known atomic masses, it is by no means obvious whether the particle number dependence 
obtained from global calculations at extreme values of N/Z is correct. Apart from strong 
theoretical and experimental interest in nuclear physics aspects of exotic nuclei, calculations 
for nuclei far from stability have strong astrophysical implications, especially in the context 
of the r-process mechanism |jT2|,|l3|. 

In previous work several aspects of nuclear structure at the limits of extreme isospin 
were discussed by means of the macroscopic-microscopic approach. In the present study, 
ground-state properties of drip-line systems and the sensitivity of predictions to effective 
forces are investigated by means of the self-consistent Hartree-Fock-Bogoliubov (HFB) ap- 
proach. 

The paper is organized as follows. Section ^ discusses the effective interactions employed 
in this study. Since pairing correlations are crucial for the behavior of drip-line systems, par- 
ticular attention is paid to the particle-particle (p-p, pairing) component of the interaction. 
After a short review of general properties of effective pairing interactions, with emphasis on 
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the density dependence, the pairing forces investigated in our work, namely contact forces 
(delta interaction, density-dependent delta interaction, and Skyrme interaction) and the 
finite-range Gogny force, are described. 

The basic ingredients of the HFB formalism in the coordinate representation (single- 
quasiparticle orbitals, time-reversal symmetry, canonical states, and various densities) are 
defined in Sec. |T|. In contrast to the single-quasiparticle wave functions which often contain 
a scattering (outgoing) component, canonical states (Sec. [Ill B\) are always localized, even 
if they have positive average energy. The interpretation of particle and (especially) pair 
densities in terms of single-particle and correlation probabilities is given in Sec. [Ill (J| . This 
interpretation is essential when relating the calculated HFB densities and fields to various 
experimental observables. 

The structure of the HFB equations is analyzed in Sec. |IV|. Here, various functions enter- 
ing the equations of motion (i.e., mass parameters and mean field potentials) are introduced 
for both p-h and p-p channels (Sees. |1V A| and |1V B| ). 

The advantage of using the coordinate-space HFB formalism for weakly bound systems 
is that in this method the particle continuum is treated properly. This important point is 
discussed in detail in Sec. 0. In particular, the difference between the single-particle Hartree- 
Fock (HF) spectra and canonical HFB spectra (Sec. |V B| ), the asymptotic properties of the 



HFB states (Sec. |V E|) and densities (Sec. |V b] } , and the effect of the pairing coupling to 
positive-energy states (Sec. |VG| ) are carefully explained. 

The robust predictions of the formalism for various experimental observables (pairing 
gaps and pair transfer amplitudes, masses and separation energies, radii, shell gaps, and 
shell structure) are reviewed in Sec. |VI|, where experimental fingerprints of the surface- 
peaked pairing fields and the quenching of shell effects far from stability are also given. 
Section |V11| contains the main conclusions of the paper. The technical details (i.e., the form 
of a mean-field Gogny Hamiltonian and the discussion of the energy cut-off in the Skyrme 
model) are collected in the Appendices. 



II. EFFECTIVE INTERACTIONS IN THE P-P CHANNEL 

The uniqueness of drip-line nuclei for studies of effective interactions is due to the very 
special role played by the pairing force. This is seen from approximate HFB relations 



between the Fermi level. A, pairing gap. A, and the particle separation energy, S [15| 



5~_A-A. (2.1) 

Since for drip-line nuclei S is very small, A + A^O. Consequently, the single-particle field 
characterized by A (determined by the p-h component of the effective interaction) and the 
pairing field A (determined by the p-p part of the effective interaction) are equally im- 
portant. In other words, contrary to the situation encountered close to the line of beta 
stability, the pairing component can no longer be treated as a residual interaction; i.e., a 
small perturbation important only in the neighborhood of the Fermi surface. 

Surprisingly, rather little is known about the basic properties of the p-p force. In most 
calculations, the pairing Hamiltonian has been approximated by the state-independent se- 
niority pairing force, or schematic multipole pairing interaction [|TB[. Such oversimplified 
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forces, usually treated by means of the BCS approximation, perform remarkably well when 
applied to nuclei in the neighborhood of the stability valley (where, as pointed out above, 
pairing can be considered as a small correction). As a result, considerable effort was devoted 
in the past to optimizing the p-h part of the interaction, while leaving the p-p component 
aside. 

Up to now, the microscopic theory of the pairing interaction has only seldom been ap- 
plied in realistic calculations for finite nuclei (see Ref. WTl for a recent example). A "first- 



principle" derivation of pairing interaction from the bare A^A^ force using the renormalization 
procedure (G- matrix technique), still encounters many problems such as, e.g., treatment of 
core polarization [|T8],|l9|. Hence, phenomenological pairing interactions are usually intro- 
duced. Two important open questions asked in this context are: (i) the role of finite range, 
and (ii) the importance of density dependence. Since the realistic effective interactions are 
believed to have a finite range, the first question seems purely academic. However, the re- 
markable success of zero-range Skyrme forces suggests that, in many cases, the finite-range 
effect can be mocked up by an explicit velocity dependence. To what extent this is true for 
the pairing channel remains to be seen. One obvious advantage of using finite-range forces 
is the automatic cut-off of high-momentum components; for the zero-range forces this is 
solved by restricting the pair scattering to a limited energy range and by an appropriate 
renormalization of the pairing coupling constant (see Appendix 0) . 

The answer to the question on the density dependence is much less clear. Early calcu- 
lations pO|,pT| for nuclear matter predicted a very weak ^5*0 pairing at the saturation point 



(/ci?=1.35fm^^). Consequently, it was concluded that strong pairing correlations in finite 
nuclei had to be due to interactions at the nuclear surface. This led to the surface delta 
interaction (SDI) |Q, a highly successful residual interaction between valence nucleons. 
Of course, the SDI is an extreme example of surface interaction. More realistic density- 
dependent pairing forces are variants of the density-dependent delta interaction (DDDI) 
introduced in the Migdal theory of finite Fermi systems . 

Since the effective interactions commonly used in the HF calculations are bound to be 



density dependent in order to reproduce the compressibility of the infinite nuclear matter ||24 
(an explicit density dependence is also said to account for three- and higher-body components 
of the interaction), it seems natural to introduce the density dependence in the p-p channel 
as well |2§. 

Interestingly, the presence (absence) of the density dependence in the pairing channel 
has consequences for the spatial properties of pairing densities and fields. As early recog- 



nized [26 1, the density- independent p-p force gives rise to a pairing field that has a volume 



character. For instance, the commonly used contact delta interaction, 

V\r,r') = Vo6{r-r'), (2.2) 
leads to volume pairing. By adding a density- dependent component, the pairing field be- 



comes surface-peaked. A simple modification of force (|2.2| ) is the DDDI P5| , |27| , |28 



V'P{r, r') = Vo6{r - r') {1 - [p(r)/pj^} , (2.3) 

where p(r) is the isoscalar nucleonic density, and Vo,Pc and 7 are constants. If pc is chosen 
such that it is close to the saturation density, pc^p{r = 0), both the resulting pair density 
and the pairing potential A(r) (see Sees. pTCl and |iVA| ) are small in the nuclear interior. 
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By varying the magnitude of the density-dependent term, the transition from volume pairing 
[pc^p(O)] to surface pairing can be probed. 

What are the experimental arguments in favor of surface pairing? Probably the strongest 
evidence is the odd-even staggering in differential radii, explained in terms of the direct cou- 
pling between the proton density and the neutron pairing tensor p9|-p^. Other experimental 
observables which strongly reflect the spatial character of pairing are the particle widths and 
energies of deep- hole states PB| , |5^ , strongly influenced by the pairing- induced coupling to 
the particle continuum, and the pair transfer form factors, directly reflecting the shape of 
the pair density. Because of strong surface effects, the properties of weakly bound nuclei are 
sensitive to the density dependence of pairing. In particular, the same type of force is used 
to describe the spatial extension of loosely bound light systems p5H38|l. (The measurable 
fingerprints of surface pairing in neutron-rich systems are further discussed in Sec. |VT[) In 
this context, it is also worth mentioning that the self-consistent model with the DDDI has 
recently been used to describe the nuclear charge radii and the moments of inertia of 
superdeformed nuclei EO]. In the latter case, the inclusion of a density dependence in the 



p-p channel turned out to be crucial for the reproduction of experimental data around ^^^Hg. 

In a series of papers |]4T|- ^] the Quasiparticle Lagrangian method (QLM) ||4^ based on 
the single-particle Green function approach in the coordinate representation has been 
applied to the description of nuclear superfluidity. The resulting pairing interaction, based 
on the Landau- Migdal ansatz [^, has zero-range and contains two-body and three-body 
components, thus leading to a density- dependent contact force similar to that of Eq. (^.31). 
(Note that in the approximations of Ref. I^TJ, the neutron pairing interaction is proportional 
to the proton density and vice versa.) However, in practical QLM calculations a 
pure density-independent delta force was used. 

A better understanding of the density dependence of the nuclear pairing interaction is 
important for theories of superfluidity in neutron stars. As pointed out in Ref. [^], it is 
impossible at present to deduce the magnitude of the pairing gaps in neutron stars with 
sufficient accuracy. Indeed, calculations of ^5*0 pairing gaps in pure neutron matter, or 
symmetric nuclear matter, based on bare NN interactions |47| suggest a strong dependence 
on the force used; in general, the singlet-S* pairing is very small at the saturation point. On 
the other hand, nuclear matter calculations with an effective finite-range interaction, namely 
the Gogny force [^], yield rather large values of the pairing gap at saturation (A~0.7MeV). 
(For relativistic HFB calculations for symmetric nuclear matter, see Ref. [|l8l. The pairing 
properties of the Skyrme force in nuclear matter were investigated in Ref. 



JSI . See also Ref. 



|50| for schematic calculations of pairing properties in nuclear matter based on the Green 
function method with a contact interaction, and Ref. |^ for a semi-classical description of 
neutron superfiuidity in neutron stars using the Gogny force.) 

In this study, several self-consistent models based upon the HFB approaches are used. 
The effective interactions employed, and other model parameters, are briefly discussed below. 

The spherical HFB-Skyrme calculations have been carried out in spatial coordinates 
following the method introduced in Ref. |^ and discussed in detail in Sees. fT|-|V]. Several 
effective Skyrme interactions are investigated. These are: (i) the Skyrme parametrization 
SkP introduced in Ref. |Q (SkP has exactly the same form in the particle-hole (p-h) and 
pairing channels); (ii) Skyrme interaction SkP*^ of Ref. (in the p-h channel, this force is 



the SkP Skyrme parametrization, while its pairing component is given by delta interaction, 
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Eq. ( p.2|) ; (iii) the Skyrme interaction SkP'^'' of Ref. |^ (in the p-h channel, this force is 
the SkP Skyrme parametrization, while its pairing component is given by Eq. (|2.3| ); (iv) 
the force SIII"^ (in the p-h channel, this is the SIII Skyrme parametrization its pairing 
component is given by the delta force of Ref. [^); (v) the force SkM*^ (in the p-h channel, 
this is the SkM* Skyrme parametrization , and its pairing part is given by the delta force 
with the parameters of Ref. 



Apart from other parameters, the above Skyrme forces differ in their values of the effective 
mass for symmetric nuclear matter, m*/m. Namely, m* /m is 0.76, 0.79, and 1 for SIII, SkM*, 
and SkP, respectively. All HFB-Skyrme results have been obtained using the pairing phase 
space as determined in Ref. |Q (see also discussion in Appendix |B|). 

A set of spherical HFB calculations has also been performed using the finite-range 
density-dependent Gogny interaction DIS of Ref. |]56|. In this effective interaction 
the central part consists of four terms parametrized with finite-range Gaussians (see Ap- 
pendix ^ . Spin-orbit and density-dependent terms of zero range are also included as in the 
Skyrme parametrizations. The pairing field is calculated from the DIS force, i.e., the same 
interaction is used for a microscopic description for both the mean field and the pairing 
channels. However, by a specific choice of the exchange contribution, the pairing component 
of the DIS is density- independent. It is also interesting to note that the pairing component 
of the DIS is repulsive at short distances and attractive at long ranges [18,58|. For the DIS 
force, the effective mass for infinite nuclear matter is m*/m=0.70. 

The parameters of the DIS interaction were chosen to reproduce certain global properties 
of a set of spherical nuclei and of nuclear matter [^. The HFB+Gogny results presented 



here were obtained by expanding the HFB wave functions in a harmonic oscillator basis 
containing up to 19 shells. 



III. INDEPENDENT-QUASIPARTICLE STATES 

The HFB approach is a variational method which uses independent-quasiparticle states 
as trial wave functions. These states are particularly convenient when used in a variational 
theory, because, due to the Wick theorem one can easily calculate for them the average 
values of an arbitrary many-body Hamiltonian. Even if the exact eigenstates of such a 
Hamiltonian can be rather remote from any one of the independent-quasiparticle states, one 
can argue ||61| that one may obtain in this way fair estimates of at least one-body observables. 



An independent-quasiparticle state is defined as a vacuum of quasiparticle operators 
which are linear combinations of particle creation and annihilation operators. This linear 



combination is called the Bogoliubov transformation ||62|- |6^] . According to the Thouless 
theorem [^], every independent-quasiparticle state |\I'), which is not orthogonal to the 
vacuum state |0), i.e., (0|\I')7^0, can be presented in the form 

l^) = exp|-^E^>X}|0), (3.1) 

where the Thouless matrix Z is antisymmetric, = —Z*, and in general complex. The 
phase of state ( p.l| ) is fixed by the condition (0|\1/)=1; the norm is given by (\l/|\l/)=det(l + 
Z'^Zy^'^. In the following, state will represent the I^=0~^ ground state of the even-even 
system. 
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We refer to standard textbooks ||2^ for a discussion of the properties of the Bogohubov 
transformation. Here we start our discussion from the trial wave function (|3.1| ) which is 
parametrized by the matrix elements of Z . This form of the independent-quasiparticle state 
is very convenient in variational applications because variations with respect to all matrix 
elements Z^^v=—Zy^ are independent of one another. 

Instead of using the matrix representation corresponding to a set of single-particle cre- 
ation operators a+ numbered by the discrete index /x, one may use the spatial coordinate 
representation. This is particularly useful when discussing spatial properties of the varia- 
tional wave functions and the coupling to the particle continuum. Therefore, in the follow- 
ing, we shall consider the operators creating a particle in the space point r and having the 
projection of spin cr=±i, 

at = E^;M<, (3.2) 

where ^/'^(ro") is the wave function of the /x-th single-particle state. To simplify the following 
expressions, we consider only one type of particle. A generalization to systems described by 
a product of neutron and proton wave functions is straightforward, while that involving the 
mixing in the isospin degree of freedom is discussed in Ref. |6^ . 
The inverse relation with respect to ( |3.2| ) is given by 



[ d'rJ2MrcT)4.- (3.3) 



Equations ( p.2[ ) and (|3.3|) assume that the wave functions ip^{ra) form an orthonormal and 
complete set. In practical calculations, the basis has to be truncated and the completeness is 
realized only approximately. The choice of the single-particle wave functions used (size of the 
set and, in particular, the asymptotic behavior) is of crucial importance to the phenomena 
discussed in this study. 

In coordinate space, the Thouless state ( p.l| ) has the form 

V), (3.4) 



1^) = exp I / dVdV J2 Z^{ra, r'a') 



and is defined by the antisymmetric complex function, Z+(ro", r'cr')= — Z*(r(T, rV), of 
space-spin coordinates. Already, at this point, we see that any variational method em- 
ploying an attractive effective interaction for a hound finite system must lead to functions 
which are localized in space, 

lim Z{ra^ r'a') = 0, for any r', a', and a. (3.5) 

|r|^oo 

Recall that in the coordinate space, values of the function Z{ra, r'a') at different space-spin 
points are the variational parameters, and that any arbitrarily small value of this function at 
large distance, \r\ — >■ cxo, would create at this point a non-zero probability density. Whether 
this would be energetically favorable depends upon the number of particles in the system 
and on the interaction used in the variational method. Apart from exotic phenomena such as 
halos, and apart from infinite matter such as in the neutron-star crust, we assume that the 
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attractiveness of the interaction always favors compact, localized probability densities, and 
hence we require the localization condition (|3.5| ) for the variational parameters Z{ra,r'a'). 

An expansion of the variational function Z{r(T, r'a') in terms of the single-particle wave 
functions is a straightforward consequence of transformations ( p.2|) and (|3.3| ), 

Zira,r'cr') = J2r,ira)Z,,r.ir'a'). (3.6) 

The localization condition, Eq. (|3.5|) , can, therefore, be guaranteed in the most economic 
way by requiring that all single-particle wave functions ip^{ra) vanish at large distances. 
Of course, this is only a matter of convenience and manageability, because any localized 
function can be expanded in any complete basis. It is, however, obvious that such an 
expansion converges very slowly if the basis has inappropriate asymptotic properties. For 
example, one can expect that a plane- wave expansion of Z{ra, r'a') would require an infinite 
number of basis states ip^lra), and in practice, any reduction to a finite basis would lead 
to serious errors. A discussion pertaining to asymptotic properties of functions in spatial 
coordinates, and the choice of an appropriate single-particle basis, will be a pivotal point in 
our study. 



A. Time-reversal 



The present study is entirely restricted to an analysis of ground-state phenomena, and 
therefore, we use only time-even variational independent-quasiparticle wave functions. The 
time-reversal operator can be represented as a product of the spin-flip operator and the 
complex conjugation; 



I.e., 



T-- 



-icTyK 



67|. The explicit time-reversed creation operators 



then have the form 



(3.7a) 
(3.7b) 



We now suppose that the set of basis states represented by the creation operators is 
closed with respect to time reversal, and that the state T+a+T is actually proportional (up 
to a phase factor Sp,=—Sfj,, |s^|=l) to another basis state denoted by a bar over the Greek 
symbol, i.e., 

T^a'j;f = spa-^, (3.8a) 
Sf,ij-,{ra) = 2ar,{r,-a). (3.8b) 

In this way, the single-particle basis is assumed to be composed of pairs of time-reversed 
states denoted by indices fj, and /2. In what follows, we use the convention that /i=/i, and 
that the sums over either /i or /i are always performed over all basis states. The phase factors 
depend on relative phases chosen for the /i-th and fi-th states of the basis; it is convenient 
to keep them unspecified in all theoretical formulae and to make a definite suitable choice 
of the phase convention only in a specific final application. 
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B. Canonical basis 



A requirement of the time- reversal symmetry of the quasiparticle vacuum or ( ^.4[ ), 
T'|\l/) = |\l/), leads to the following conditions: 

= s;s:z;„ (3.9a) 
Z{ra, r'a') = Aaa'Z*{r, -a, r', -a'). (3.9b) 

These properties allow the introduction of more suitable forms of Z^j^ and Z{ra,r'a'); 
namely, 

Zfj,u '■= Sfj,Zfiiy, (3.10a) 

Z{ra, r'a') := 2aZ{r, -a, r'a'). (3.10b) 

The matrix Z^^^ and the function Z{ra, r'a') are both time-even and hermitian, 

= (3.11a) 
Z*{ra, r'a') = Z{r'a', ra), (3.11b) 

and therefore they can be considered as usual operators in the corresponding Hilbert spaces. 
In particular, the function Z*{ra, r'a') can be diagonalized by solving the following integral 
eigenequation: 

/ dV^Z(ra,rV)4(rV) = z^i>^{ra), (3.12) 

a' 

where are real eigenvalues, z^=Zfi. The eigenfunctions ip^{ra) form the single-particle 
basis, usually referred to as the canonical basis. Canonical states, together with the eigen- 
values Zfj^, completely define the quasiparticle vacuum |\E'). (Here and in the following we use 
the checked symbols, e.g., ip^ and a^, to denote objects pertaining to the canonical basis.) 

Two important remarks concerning the canonical basis are now in order. First, the 
localization condition ( |3.5| ) directly results in the fact that all canonical-basis single-particle 
wave functions ip^{ra) are localized in space; i.e., vanish at large distances, \r\ oo. 
Therefore, as discussed previously, a choice of the localized wave functions for the basis 
states il)^j_{ra) may allow for a rapid convergence in the expansion 

^^(ra) = ^Dy^,i)^{ra), (3.13a) 

a'^ = J2Dui,at- (3.13b) 

V 

Second, since Z{ra, r'a') and Z^y are related by [cf. Eq. (|3.6|) ] 

~Z{raya') = Y.i,^{ra)Z^,i,l{r' o'). (3.14) 

a diagonalization of Z{ra, r'a'), Eq. (|3.12|) , is equivalent to a diagonalization of Z^^, 
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'^^Z^yDyr — ZrD^T-. (3.15) 



Therefore, in the canonical basis, the Thouless state ( |3.1| ) acquires the well-known separable 
BCS-like form 

1^) = exp|i^Z^^s^ata+| |0) 



exp i ^ , 

ly>0 



z^Syc4a^ I |0) 

\{{l + z,sXK)\^). (3.16) 



u>0 



where the symbol > denotes the sum over one-half of the basis states with only one 
state (either one) of each time-reversed pair (z/, v) included, and is the creation operator 
in the canonical basis. 



C. Density matrices and the correlation probability 

According to the Wick theorem [|60ip4| for the independent-quasiparticle state, Eqs. 
( p.l|) or (|3^), an average value of any operator can be expressed through average values of 
bifermion operators, 

p{ra,r'a') = (^la^'flral^), (3.17a) 
p{ra,r'a') = -2(T'(^|a^/,_<,/a^<,|^). (3.17b) 

Functions p{ra, r'a') and p{ra, r'a') are called the particle and pairing density matrices, 
respectively. For a time-reversal invariant state |\E'), both density matrices are time-even 
and hermitian: 

p{ra, r'a') = Aaa' p{r -a,r' - a')*, (3.18a) 
p{ra, r'a') = Aaa'p{r -a,r' - a')*. (3.18b) 

Therefore, the pairing density matrix p{ra, r'a') is more convenient to use than the standard 
pairing tensor K,{ra,r'a') 

K{ra, r'a') = 2a'p{ra, r', -a'), (3.19) 

which is an antisymmetric function of the space-spin arguments. 

The formulae expressing p{ra, r'a') and p{ra, r'a') in terms of the function Z{ra, r'a') 
can be easily derived from those for the density matrix and the pairing tensor , and they 
read 

p={l + ZY^Z\ (3.20a) 
p={l + Z'^)-^Z. (3.20b) 

As a result, the density matrices obey the following relations: 
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P- P = P- P, (3.21a) 
p ■ p + p ■ p = p. (3.21b) 

In the above equations, the matrix multiplications and inversions should be understood in 
the operator sense; i.e., they involve the integration over space and summation over spin 
variables. For instance: 

{p- p){ra,r'a') = f dV J] p(r(T, rV')p(r'V', rV). (3.22) 

Local HFB densities, i.e., the density matrices for equal spatial arguments, r'=r, have 
very well-defined physical interpretations. To see this, let us assume that ipxs{rcr) is a 
normalized single-particle wave function (wave packet) concentrated in a small volume Vx 
around the point r=x and having the spin s=a. The corresponding creation operator 

<s= [ dV^^,,(ra)a+ , (3.23) 

together with its hermitian conjugate, define the operator 

Nxs = a+ (3.24) 

which measures the number of particles in the vicinity of the point x. Since 

NL = Nxs, (3.25) 

Nxs can be regarded as a projection operator which projects out the component of the 
many-body wave function that contains one spin-s fermion in the volume Vx- Therefore, its 
average value gives the probability to find a particle with spin s in this volume: 

Viixs) = i^lNxsl"^) = Vxp{xs, xs). (3.26) 

In a very similar way, the probability of finding a fermion in Vx having opposite spin can 
be obtained by considering the time-reversed wave function 2cr'?/'*g(r, — cr), cf. Eq. ( |3.7bD . 
This gives 

Vi{x, -s) = (^|iV,,_,|^) = Vxp{x, -s, X, -s). (3.27) 

Due to time-reversal symmetry, probabilities ( |3.26|) and ( p. 271) are equal. 

We may now ask the question, "What is the probability of finding a pair of fermions 
with opposite spin projections in the volume Vx, V2{x)T\ If one considers two indepen- 
dent measurements, where in the first one is found the spin-s fermion, and in another 
one the spin-(— s) fermion, V2{x) is equal to the product of individual probabilities; i.e., 
Vi{x, s)Vi{x, —s). On the other hand, if one wants to find in Vx both fermions simulta- 
neously, one should project out from a corresponding two-fermion component. In this 
case, V2{x) becomes the expectation value of the product of the projection operators A^a;^ 
and Nx-s', i-e., V2{x) = {"^iNxsNx-sl"^)- Using the Wick theorem, this average value is 

'P2ix) = V^p{xs, xs)p{x, -s, X, -s)+ ,g 
Vxpixs,xs)p{x,-s,x,-s), 
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or in terms of the time-even spin-averaged densities: 



for 



V2{x) = \vlp{xf + \v^~p{x)\ (3.29) 



P(^) = H p{r(y, ra), (3.30a) 

cr 

p{r) ='^p{ra,ra). (3.30b) 



Since the first terms in Eqs. ( |3.28D and (|3.29|) describe the probabihty of finding the two 
fermions in independent measurements, the second terms in these equations should be in- 
terpreted as the probabihty of finding the correlated pair at point x. 

The above arguments aUow us to give a transparent physical interpretation to the local 
HFB densities. Namely, as usual, p{r) represents the probability density of finding a particle 
at the given point. On the other hand, p(r)^ gives the correlation probability density; i.e., 
the probability of finding a pair of fermions in excess of the probability of finding two 
uncorrelated fermions. 

It is important to note that kinematic conditions (|3.21|) , which result from the fact that 
|\E') is an independent-quasiparticle state, Eq. ( p.4|) , do not directly constrain the local values 
of the particle and pairing density matrices. In particular, there is no obvious kinematic 
relation between the probability of finding two independent particles at a given point of 
space, and the probability of finding a correlated pair at the same point. In particular, the 
first one can be small, while the second can be large (see discussion in Sees. [Ill C 1 and 



|V b] ) . This result means that in such a situation the experiments probing the presence of 
two particles will always find these two particles as correlated pairs without a "background" 
characteristic of two independent particles. 

Relations ( |3.20| ) imply that all three functions: Z{ra,r'a'), p{ra,r'a'), and p{ra,r'a') 
are diagonal in the canonical basis, cf. Eq. ( p.l2| ). Using the standard notation for the 
eigenvalues of p and p, one obtains 

/ d'r'J2pira,r'a')i>,ir'a') = vl^,{ra), (3.31a) 

a' 

I dV'5]p(ra,rV)V^^(rV) = u^v^^^{ra), (3.31b) 

a' 

where the real factors and are given by 

A completeness of the canonical basis leads to standard expressions for the density matrices: 

p(ra,rV) = J] t;JV^;(ra)^^(rV), (3.33a) 
p{ra,r'a') = ^M^f^^*(ra)V^^(rV'), (3.33b) 
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Equation ( p.31a|) represents the traditional definition of tlie canonical states as the eigen- 
states of the HFB density matrix. It also shows that the canonical states are the natural 



states [|68|-[r3[1 for the density matrix corresponding to the independent-quasiparticle many- 
body state Eq. (|3.17a|) , and the eigenvalues are the corresponding natural occupation 
numbers. 

One may now easily repeat the previous analysis of probabilities of finding a particle, or 
a pair of particles, in the canonical-basis single-particle state 'ip^{ra). The result, analogous 
to Eqs. ( p. 2^ ) and (|3.29| ), is Vi{fi)=v^, and 'P2(/^)=w^f^+f^. In this case, due to the 
normalization condition u^+f^=l, Pi(/i)='P2(/^)- This result means that the particles in the 
canonical states with indices /i and ft are extremely correlated spatially; i.e., the probability 
of finding the canonical pair, w^f^, is directly dependent on the probability of finding two 
independent canonical fermions, fj^. However, as discussed above, a similar direct relation 
between Vi{x) and V2{x) does not exist. In particular Vi{x)^V2{x). 



Examples of particle and pairing densities 



Figures and display the particle and pairing local spherical neutron HFB densities 
p(r) and p(r), Eq. ( 3.3(J| ), as functions of the radial coordinate r=|r|. Results are shown for 



several tin isotopes across the stability valley. For particle densities, the results obtained 
with the SkP and SkP*^ interactions are almost indistinguishable. Therefore, Fig. (middle 
panel) shows results for the SIII*^ interaction. For pairing densities, compared in Fig. ^ are 
results for SkP, SkP*^, and DIS effective interactions. 

The particle densities obtained with these three effective interactions are qualitatively 
very similar. One can see that adding neutrons results in a simultaneous increase of the 
central neutron density, and of the density in the surface region. The relative magnitude 
of the two effects is governed by a balance between the volume and the surface asymmetry 
energies of effective interactions. Since all three forces considered have been fitted in a similar 
way to bulk nuclear properties, including the isospin dependence, the resulting balance 
between the volume and the surface isospin effects is similar. Of course, this does not 
exclude some differences which are seen when a more detailed comparison is carried out. 

The pairing densities shown in Fig. ^ reflect different characters of the interactions used 
in the p-p channel. The contact force (the SkP*^ results) leads to the pairing densities which 
are, in general, largest at the origin and decrease towards the surface. (This general trend 
is slightly modified by shell fluctuations resulting from contributions from orbitals near the 
Fermi level.) At the surface, the isospin dependence of SkP*^ is fairly weak. For example, 
there is very little difference between the pairing densities in ^^°Sn and ^''^Sn. These results 
are characteristic for the volume-type pairing correlations. 

A different pattern appears for the SkP results, where the density dependence renders 
the p-p interaction strongly peaked at the surface. In this case, the pairing densities tend 
to increase when going from the center of the nucleus towards its surface. Again, the shell 
fluctuations are superimposed on top of this general behavior. In particular, the central 
bump in the pairing density in ^^''Sn is due to a contribution from the 3si/2 state. A more 
pronounced dependence on the neutron excess is seen in the surface region. Especially near 
the drip line, the pairing density develops a long tail extending towards large distances. 
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The results obtained for the finite-range interaction DIS exhibit intermediate features 
between the surface and the volume type of pairing correlations. In particular, in the nuclear 
interior one observes a fairly large region of relatively constant pairing density. The overall 
magnitude of the pairing densities is very similar in all three approaches. In particular, it 
is interesting to see that at the nuclear surface (r~5fm) all three pairing densities in i20Sn 
are very close to 0.018 fm~^. 



IV. HARTREE-FOCK-BOGOLIUBOV EQUATIONS 



We begin this section by presenting basic definitions and equations of the HFB approach. 
The HFB theory is discussed in many textbooks and review articles (see Refs. |7^J2^, for 
example), while its aspects pertaining to the coordinate representation have been presented 
in Ref. |^2[|. An earlier discussion of the coordinate-representation HFB formalism has 
been given by Bulgac, whose work is available only in preprint form |Q. Recently, similar 
methods have also been applied to a description of light nuclei |]37|j38| . It is also worth 



mentioning that the Green function approach in the coordinate representation (the Gor'kov 
method ||75|), is formally equivalent to HFB - cf. discussion in Refs. [0,0. The only 



difference between the methods lies in the explicit energy dependence of the quasiparticle 
mass operator, an analog to the p-h single-particle HF Hamiltonian (see below). 



A. HFB energy and HFB potentials 

The two-body effective Hamiltonian of a nuclear system can be written in the coordinate 
representation as 

H = J Srd'^r'J2T{r(T,r'(T')ataar'a' (4.1) 

ctct' 

+ i I dVidV^dV^dV'^ ^2^2; r[a[, ry,)a^^^ 



„ „ dr' rr' CLf' rr' ■ 

1(71 r2(T2 ^2^2 



(Tia2Cf-^(T2 



The first term represents the kinetic energy, while the second one is the two-body interaction. 
In the following, we assume that V{riai, r2(72] '''2^2) includes the exchange terms. 



The average energy of the Hamiltonian (|4.1| ) in a time-even HFB vacuum (O) reads 
Ehfb = J d^rd^r'J2T{ra,r'a')p{r'a',ra) (4.2) 

era' 

+ \j dVidVadV'idVa ^ l^(riai, ra^s; r\o[, r'^a'^) p{r\a[, ri(Ti)p(r'2(T^, ra^s) 

~\j d^'^idVadV^dVa ^1 4(Ti(T2V(ri, -CTi, ra^a; r\o[, r'^, -a2)p(riai, r2(r2)p{r[a[, rsag)- 



ai(T2(T-^cr2 



The last two terms are the interaction energies in the particle-hole (p-h) and in the particle- 
particle (p-p) channels, respectively. Equivalently, one can define the p-h and p-p single- 
particle Hamiltonians, h{ra,r'a') = T{r(j,r'a') + T{ra,r'a') and h{ra,r'a'), respectively: 
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r(ra, r'a') = J d^r^d'r', ^ V{ra, r^a^, r'a\ r',a',)p{r',a',, rscxs), (4.3a) 
h{ra, r'a') = [ dV;dV^ ^ 2a'a',V{ra, r', -a'; r',a'„ r'„ -a',)p{r',a'„ r',a',), (4.3b) 

which gives the HFB energy in the form: 

1 /" ~ 

)p(rV, ro) + /i(ro-, r'a')p{r'a' , ra) + /i(ro-, r'a')p{r'a', ra 

era' 

(4.4) 

Additional terms coming from the density- dependence of the two-body interaction V have 
been for simphcity omitted in Eqs. ( [4.3a| ), ( [4.3b|) , and (|4.4|) . The last term in Eq. 



impair = \j dVdV hira, r'a')pir'a', ra), (4.5) 

era' 

represents the pairing energy. We also define the average magnitude of pairing correlations 
by the formula 

(A) = J dVdV J2 r'a')p{r'a', ra), (4.6) 

aa' 

where N'^ is the number of particles (neutrons or protons). 

The p-h and p-p mean fields ( [4.3[ ) have particularly simple forms for the Skyrme inter- 
action In Appendix ^ we present the form of the p-h and p-p mean-field Hamiltonians 
in the case of a local two-body finite-range Gogny interaction. 

1. Examples of the p-h and p-p potentials 

In this section we aim at comparing the self-consistent potentials obtained with the 
Skyrme and Gogny forces. Such a comparison cannot be carried out directly, because the 
corresponding integral kernels h{ra, r'a') and h{ra, r'a') have different structure. For the 
Skyrme interaction, they are proportional to 6{r—r') and depend also on the differential 
operators (linear momenta) [Q, while for the Gogny interaction they are sums of terms 
proportional to 6{i — r') and terms which are functions of r and r' (Appendix 

Therefore, for the purpose of the present comparison we introduce operational prescrip- 
tions to calculate the local parts of the integral kernels: 

U{r) = COC [r{ra, r'a')] , (4.7a) 
U{r) = COC [h{ra, r'a')] . (4.7b) 

These formal definitions in practice amount to: (i) disregarding the momentum-dependent 
terms of the kernels, (ii) considering only terms with a=a'=l/2 (which by time reversal 
symmetry are equal to those with a=a'=— 1/2), and (iii) taking into account only the term 
proportional to 6{r—r'), if such a term is present. The expressions for U{r) and U{r) 
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can be found in Appendix A of Ref. [|52| (Skyrme interaction) and in Appendix A (Gogny 
interaction). In the Skyrme calculations, the contribution of the Coulomb interaction to 
U (r) has been neglected since it is estimated to be small. 

In the case of finite-range local interactions (such as Gogny or Coulomb), the correspond- 
ing non-local pairing field h{ra, r'a') does not contain the term proportional to 6{r—r') (see 
Appendix Consequently, the local field U{r) cannot be extracted in a meaningful way. 
For instance, the diagonal (i.e., r'=r) part of of the DIS pairing field is positive; i.e., it 
is dominated by the short-range repulsive component rather than the long-range attractive 
part 

In the spherical case, the potentials U (r) and U (r) depend on only one radial coordinate 
r=|r'|. This facilitates the qualitative comparison between different forces. Figure ^displays 
the self-consistent spherical local p-h potentials U{r), Eq. ([4.7|), for several tin isotopes, 
calculated with SkP, SIII"^, and DIS interactions (the results with SkP*^ are very close to those 
with SkP). The terms depending on the angular momentum, which result from a reduction 
to the radial coordinate, are not included. (The general behavior of the self-consistent p-h 
potentials has already been discussed many times in the literature, e.g. ||76|^r8[, and we 
include these results only for completeness and for a comparison with the corresponding p-p 
potentials, for which the detailed analysis does not exist.) 

Qualitatively, the results for U{r) obtained with different effective forces are quite similar, 
which reflects the fact that all these interactions correctly describe global nuclear properties. 
In particular, one sees that with increasing neutron excess the neutron potentials become 
more shallow in the interior and more wide in the outer region. Interestingly, for each of 
these three forces there exists a pivoting point at which the potential does not depend on the 
neutron excess. For the three forces presented, this occurs at r=5.9, 4.6, and 5.4 fm, respec- 
tively. The differences in the overall depths of the average potentials reflect the associated 
effective masses (i.e., the non-local contributions of the two-body interactions). 

The analogous results for the p-p potentials f/(r) calculated for the SkP and SkP'' inter- 
actions are shown in Fig. ^. On can see that the different character of pairing interactions is 
directly reflected in the form of the p-p potentials. Particularly noteworthy is the fact that 
the density-dependent pairing interaction in SkP yields a very pronounced surface-peaked 
potential (the behavior of f/(r) at large distances is further discussed in Sec. |VG| ). One 
can easily understand its form by recalling that this potential is equal to the product of the 
pairing density p{r) [Fig. |^ and the function which roughly resembles the behavior of DDDI 
of Eq. ( |2.3| ); i.e., small in the interior and large in the outer region. Of course, values of p{r) 
and U{r) depend on each other by the fact that they both result from a self-consistent solu- 
tion of the complete HFB equation in which the p-h and p-p channels are coupled together 
(see Sec. [IV C| ) . Similar results were also obtained in Refs. |^ (in the HFB+SkP model) 
and 1^ (in the QLM) for the proton-rich rare-earth nuclei. 

Since the p-h channel provides the bulk part of the interaction energy, the particle 
densities p(r) closely follow the pattern of the p-h potentials (i.e., the density is large where 
the potential is deep). An analogous relation is only partly true for p(r) and U{r); i.e., 
even the dramatic surface character of the SkP p-p potential (Fig. 1) does not result in the 
pairing density being similarly peaked at the surface. Recall that the contributions to p(r) 
come mainly from a few wave functions near the Fermi surface, and that the form of these 
wave functions is mainly governed by the p-h channel. Since these wave functions must 
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have significant components in the interior, the resulting pairing densities cannot exactly 
fit into the surface-peaked p-p potentials. Nevertheless, a clear tendency towards surface 
localization is evident in Fig. |^. 

In the case of the pure contact interaction (SkP*^ calculations) the p-p potential is ex- 
actly proportional to the pairing density |Q with the proportionality constant Vo/2 equal 
to — SOMeVfm^ Therefore, the resulting potential is concentrated at the origin and 

increases towards the surface. (Early calculations of p-p potentials in the QLM with the 
density-independent delta interaction can be found in Ref. |^3[. The general behavior of 
U{r), denoted as A(r) therein, is very similar to our SkP-^ results.) 



B. HFB equations in the coordinate representation 

The variation of the HFB energy with respect to independent parameters Z{ra, r'a') 
leads to the HFB equation P^J5^ , 



h{ra,r'a') h{ra,r'a') \ ( 0i(E,rV)\ _( E + \ \( ME,ra) 
h{ra,r'a') -h{ra,r'a') J {M^^r'a') ) ~ [ E-\) [ME,ra) 



(4.8) 



where (f)i{E,ra) and (j)2{E,ra) are upper and lower components of the two-component 
single-quasiparticle HFB wave function, and A is the Fermi energy. 

Properties of the HFB equation in the spatial coordinates, Eq. (^]8|), have been discussed 
in Ref. [0. In particular, it has been shown that the spectrum of eigenenergies E is contin- 
uous for \E\>—X and discrete for \E\<—X. Since for E>0 and A<0 the lower components 
(j)2{E,ra) are localized functions of r, the density matrices. 



/oo 
dn{E)(f)2{E,ra)(j);{E,r'a'), (4.9a) 

/OO 
dn{E)<P2{E,ra)<Pl{E,r'a'), (4.9b) 

are always localized. 

For the case of a discretized continuum. Sec. |V7\| , the integral over the energy reduces to 
a discrete sum but one should still carefully distinguish between contributions coming 
from the discrete (i?„<— A) and discretized {En>—\) states. The orthogonality relation for 
the single-quasiparticle HFB wave functions reads 

/ dV ^ [0*(E„, r(T)0i(E„,, ra) + 0;(E„, ra)02(^n', ra)] = 6^^. (4.10) 

a- 

It is seen from Eq. ( |4.1CI|) that the lower components are not normalized. Their norms, 

Nn= f dV^|02(K,ra)|2, (4.11) 

define the total number of particles 
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N= /dVp(r) = ^iV„. 



(4.12) 



In the HFB theory, the localization condition ( p.5|) discussed in Sec. |T| is automatically 
guaranteed for any system with negative Fermi energy A. This allows studying nuclei which 
are near the particle drip lines where the Fermi energy approaches zero through negative 
values. 

For the Skyrme interaction, the HFB equation ([4.8|) is a differential equation in spatial 
coordinates If the spherical symmetry is imposed, which is assumed in the following, 
this equation reads 



-A 


f M 


M ^ 


1- 


dr 




-My 


' dr 




r(f)i{E, r) 
r(f)2{E, r) 



E 



r(j)i{E, r) 
r(j)2{E, r) 



(4.13) 



where M and M are p-h and p-p mass parameters, respectively, and U and U are defined 
in Sec. p^V A| . Due to the spherical symmetry, Eq. ([4.13|) is solved separately for each partial 
wave (j, ^). The potentials include also the centrifugal and spin-orbit terms, and the p-h 
mass parameter M is expressed in terms of the effective mass m*; i.e., see Ref. 



1 5^ for details. 

Before discussing the properties of the HFB wave functions, we analyze the structure 
of the spherical HFB Hamiltonian of Eq. ( [4.13| ). Figure ^ shows the behavior of M(r) 
and M(r), and U{r) and U{r) (central parts only) obtained for neutrons in ^^°Sn in the 
HFB-|-SkP model. The p-h functions, M(r) and U{r), are similar to those obtained in other 
mean-field theories. M(r) has values close to h /2m~20 MeVfm^, which corresponds to the 
value of the free nucleon mass m. In the nuclear interior, this function has slightly smaller 
values, because the effective mass m* is here slightly larger than m. This effect is due to 
the non-zero isovector effective mass of the Skyrme SkP interaction; recall that for this 
interaction the nuclear-matter value of the isoscalar effective mass is m*=m. The central 
potential f/(r) has the standard depth of about 40 MeV and disappears around r=7.5 fm. 

The form of the p-p functions, M(r) and t/(r), characterizes the pairing properties of 
the system. One may note that both these functions are essentially peaked at the nuclear 
surface. In ^^°Sn they also exhibit central bumps resulting from the fact that in this nucleus 
the neutron 3si/2 orbital is located near the Fermi surface. Values of M(r) are (in the chosen 
units) an order of magnitude smaller than those of lJ{r). This should be compared with the 
results obtained for the p-h channel, where the values of M(r) are only about a factor of 2 
smaller than those of U{r). It means that, for the SkP parametrization, the kinetic term in 
the p-p channel (which simulates the finite-range effects) is relatively less important than 
the kinetic energy term in the p-h channel. 



C. Single-quasiparticle wave functions 



This section contains the discussion of HFB wave functions 
V C 1|) , canonical-basis wave functions V'^(r) (Sec. [IV C 2 




i{E,r) and (j)2{E,r) (Sec. 
and HF+BCS wave functions 



In the following, the HFB equation ( [4.13[ ) was solved in the spherical box 
of the radius i?box=20fm for the j=l/2 and i=0 (si/2) neutron states; i.e., for vanishing 
centrifugal. Coulomb, and spin-orbit potentials. The calculations were performed for ^^''Sn. 



18 



1. Examples of the single- quasiparticle wave functions 



The neutron single-quasiparticle wave functions are presented in Fig. ^. The upper 
components r(j)i{En,r), and the lower components r(l)2{En,r), are plotted in the left and 
right columns, respectively. Because a box of a finite radius was used, the particle continuum 
is discretized. The positive quasiparticle eigenenergies En are in increasing order numbered 
by the index n, and their values are tabulated in the left portion of Table |, together 
with the norms of the lower components ( [4.11| ), A''„=47r / r'^dr\(f)2{En,r)\'^. Since the lower 
components define the particle density matrix [Eq. ([4.9a|) ] the numbers {2j+l)Nn (i.e., 2Nn 
for the j=l/2 case considered) constitute contributions of a given quasiparticle state to the 
total number of neutrons (see Eq. (|4.12| )). 



Wave functions in Fig. ^, and the entries in Table |, have been ordered from the bottom 
to the top not according to the excitation-energy index n, but rather according to numbers of 
nodes of the large component. (The large component is the lower component for hole states 
and the upper component for particle states - see Fig. ^.) The lower component of the n=8 
state is large, and it has zero nodes; therefore it is plotted at the bottom of the figure. Next 
comes the n=5 state, whose lower component has one node, and the n=l state with two 
nodes. Lower components of these three states are larger than their upper components and 
they contribute almost 2 particles each to the total number of neutrons. Consequently, these 
quasiparticle states should be associated with the lsi/2, 2si/2, and 3si/2, single-hole states. 

For all other calculated Si/2 states the upper components are larger than the lower ones, 
and these states contribute small fractions to the particle number, see Table |. Consequently, 
these quasiparticle states should be associated with the Si/2 single-particle states. The 
behavior of these wave functions differs in the nuclear interior (i.e., for r < R where i?~7 fm 
is the nuclear radius) and outside (r > R). Since the wavelength of the upper component 



is roughly proportional to 1/yEn + X — U{r), the ratio of the corresponding wavelengths 
behaves as 



where U{0) is the depth of the neutron potential well. For the Si/2 neutron states in ^^°Sn the 
excitation energy, E^ + A, can be found from Table | (A=— 7.94 MeV), and [/(0)~— 45MeV 
(see Fig. |^). 

The upper component of the n=2 state has three nodes. However, for r > R the exterior 
part of the wave function corresponds to a half- wave; i.e., it represents the lowest-energy 
discretized continuum state. Since En + X is only 0.97 MeV, the wavelength in the nuclear 
interior is ~6.5 times shorter than Aout- The next two wave functions have four and five nodes 
in their upper components. Compared to the n=2 state, they exhibit shorter wavelengths 
both outside and inside the nucleus (the corresponding excitation energies are larger), and 
the ratio Aout/Ain decreases according to Eq. ( [4.14| ). 

The quasiparticle states with n=2, 3, and 4 should be associated with the 4si/2, 5si/2, 
and 6si/2 states in the particle continuum. Of course, the values of their quasiparticle 
energies strongly depend on the size of the box, because the wavelength of their exterior 
parts will increase with increasing -Rbox (is roughly proportional to -Rbox)- 
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From the above discussion, one can see that the structure of large components resembles 
very much that of the HF wave functions. Moreover, the small components are very small 
compared to the large ones; in order to plot both of them in the same scale (Fig. |^) they have 
to be multiplied by factors from 10 to 25. Only the lowest quasiparticle state (n=l), which 
corresponds to the 3si/2 state near the Fermi surface, has the two components of a similar 
magnitude. It is to be noted, however, that the detailed structure of small components is 
decisive for a description of the pairing correlations. Indeed, both components are coupled 
in the HFB equations by the pairing fields h{ra, r'a') or U. 

In agreement with general asymptotic properties of the upper and lower components 
p3| , |52[| , one sees in Fig. ^that the lower components vanish at large distances for all quasi- 
particle states, regardless of the excitation energy. Consequently, the resulting density ma- 
trix is localized. It is interesting to observe (Table |) that the norms of the lower components 
Nn do not behave monotonically with quasiparticle energy. Namely, Nn is about 0.0002 for 
n=2; then it increases to 0.0019 at n=6, and only then it decreases to about 0.0001 at 
n=ll. This means that the pairing correlations couple states with very high quasiparticle 
excitations and short- wavelength upper components; i.e., located high up in the particle 
continuum. In the considered example, only by going to the energy region of as high as 
50 MeV is the pairing coupling to the continuum states exhausted. 

Apart from the n=l state which has the quasiparticle energy E smaller than —A, for all 
other quasiparticle states the upper components oscillate at large distances; i.e., these states 
belong to the HFB continuum. This seems natural for the 4si/2, 5si/2, and 6si/2 states 
discussed above, but it also holds for the deep-hole states lsi/2 and 2si/2. This illustrates 
the physical property of the deep-hole states that once such a state is excited, it is coupled 
to the particle continuum and acquires some particle width. Of course, before such a hole 
is created (e.g., one-quasiparticle excitation in the neighboring nucleus) the nucleus (i.e., 
quasiparticle vacuum) is perfectly particle-bound and the contributions from the deep-hole- 
like quasiparticle states to the density matrix are localized in space. 



2. Examples of the canonical-basis wave functions 

By solving the integral eigenequation for the density matrix (|3.31a|) , one obtains the 
canonical-basis wave functions tpnir). Actually, when the HFB equation ( [4.13| ) is solved by 
a discretization method on a spatial mesh, as is done here, the density matrix is represented 
by a matrix and the integral eigenequation becomes the usual matrix eigenproblem. In 
the present application to ^^''Sn, the mesh of equally spaced points with Ar=0.25fm was 
used and then the canonical-basis wave functions were obtained on the same mesh of points. 
These wave functions are plotted in Fig. ^, while other characteristics of the canonical states 
are listed on the right-hand side portion of Table |. Here the states are ordered from bottom 
to top according to their occupation probabilities vf^. 

When /i increases from 1 to 5, the number of nodes of the canonical-basis wave functions 
increases from zero to four. Therefore, these states represent the lsi/2 to 5si/2 single-particle 
states. The first three of them have large occupation probabilities f^, negative average values 

of the p-h Hamiltonian, and positive pairing gaps [see Eq. ( [4.221 )]. These states have 



all the characteristics of bound single-particle states, and their wave functions strongly 
resemble the large components of the n=8, 5, and 1 quasiparticle states shown in Fig. 0. 
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It is interesting to note that the two states /i=4 and 5 follow exactly the same pattern of 
localized wave functions, despite the positive values of e^. Therefore, these two states can 
be understood as the representatives of the positive-energy spectrum in the ground-state 
of ^^°Sn. We purposely avoid using the term "particle continuum", because these orbitals 
represent discrete and localized eigenstates of the density matrix. 

Table | shows that the occupation probabilities of the canonical-basis states with 
/x=4,. . .,7 decrease very rapidly. In fact only states with /i=4 and 5 have tangible oc- 
cupation probabilities; one can say that the remaining orbitals are entirely empty. This 
feature has to be compared with the sequence of norms of the lower HFB components, Nn, 
which do not fall down to zero at even a nearly similar pace. This demonstrates that even if 
the convergence of the HFB eigenproblem requires high quasiparticle energies, the number 
of physically important single-particle states is very restrained. Unfortunately, as discussed 
below in Sec. [IV D| , one cannot obtain the canonical-basis states without actually solving 
the HFB equations up to high energies. For /i=6 and higher, the occupation probabilities 
are so small that the numerical procedure used to diagonalize the density matrix returns 
accidental mixtures of almost degenerate eigenfunctions. This is seen in Fig. |^, where the 
wave function with fi=Q has six nodes instead of five, expected from the regular sequence. 
Also the energies are for these nearly empty states randomly scattered between 40 and 
70 MeV, while the pairing gaps are scattered around zero. 



3. Examples of the BCS quasiparticle wave functions 

The BCS quasiparticle wave functions can be obtained by enforcing the BCS approxima- 
tion on the HFB equations. This is done by setting the pairing Hamiltonian hto a, constant; 
i.e., by using M(r)=0 and ?7(r)=— 1.232 MeV. This value of U is equal to minus the HFB 
average neutron pairing gap, as defined in Eq. ( [4.6|) . As seen in Fig. p, the pattern of large 
components follows closely that obtained in the HFB method, while the shapes of small 
components are entirely different. Indeed, since in the BCS approximation lower and upper 
components are simply proportional, small and large components have the same asymptotic 
properties. This leads to serious inconsistencies, because the small lower components are 
not localized any more, and introduce an unphysical particle gas in the density matrix, while 
the small upper components are localized and the corresponding deep-hole states have no 
particle width. 



D. HFB equations in the canonical basis 



It is seen in Eqs. (12) and ( [4.3|) that the two-body interaction enters the p-h and p-p 



channels in a different way. This is particularly conspicuous when the canonical basis (|3.31|) 
is used; i.e., 

^HFB = ^f^f^^l + ^l^'^K^l 

^G^yU^V^UyVy, (4.15) 



1 

4 
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where 

= \ (W + ' (4- 16a) 

G^^u = -s*^s^V^p^p. (4.16b) 
The two-body matrix elements in the canonical basis are defined as usual: 

W^''= /dVidVadVidVa ^ V(riai, ^2^2; r^ai, r2a2)^;(ri(Ti)^*(r-2a2)^^'(r;a;)V^^'(r 

(4.17) 

Since we include in V{riai,r2(J2',r[a[,r2a2) the exchange term, the matrix V^u^iu' is anti- 
symmetric in /ii/ and in /xV. Due to the hermiticity and the time-reversal symmetry of the 
interaction, matrices F^iy and G^i/ obey the following symmetry relations, 

F=F*=F=F- = F- (4 18a) 

G^u = C*!/ = Gufi = G^u = Gfiu- (4.18b) 

The matrix F is defined by different matrix elements of the interaction than the matrix G. 
Namely, the matrix element F^j^^ represents a "diagonal" scattering of pairs of states fiu^/iu 
(or /iz/^/iz/). This type of scattering concerns a// pairs of states. The resulting contributions 



to the energy, Eq. (|4.15|) , involve the occupation probabilities of the single-particle states 
constituting each pair. On the other hand, the matrix elements G^^, represent a "non- 
diagonal" scattering of pairs of time-reversed states uu^/ijl. This scattering concerns only 
a very special subset of all pairs. 

In principle, an effective interaction should describe both channels of interaction at the 
same time. This is, for example, the case for the Gogny interaction |^ and for the Skyrme 



SkP interaction ||52|. However, the fact that both channels of interaction play a different 
role in the HFB theory allows the use of different forms of interaction to model the p-h and 
p-p channels. Such an approach is additionally motivated by the fact that the interaction in 
the p-h channel, which defines, e.g., the saturation properties, is much better known than 
the p-p interaction. Moreover, the p-h channel provides a two-orders-of-magnitude larger 
interaction energy. 

Since the canonical-basis wave functions ■?/'(rcr) are all localized, it is instructive to con- 
sider the HFB equations in this particular basis. They read: 

{h - A)^^r/^^ + h^^^^^ = 0, (4.19a) 

{h-\)^,ui^,u-hf,yri^y = E^,^, (4.19b) 

where 

^n^ll' '■= u^Vu + UyV^, (4.20a) 

it,u ■■= Uf^Uu - v„Vf,. (4.20b) 



Equation ( |4.19a|) is equivalent to the variational condition that the HFB energy is mini- 



mized, while Eq. ( [4.19b|) defines the energy matrix E^^. (The matrix E^^ represents the 
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HFB Hamiltonian in the canonical basis.) Since for every pair of indices /iz/ it holds that 
(^Ati/)^+(^Mi')^~-'-' -^1^- ( l^-19|) can be written as 



-h 



(4.21a) 
(4.21b) 



The occupation probabilities are solely determined by the diagonal matrix elements of 
the p-h and p-p Hamiltonians, 



h 



and the result is 



sign(A^ 



MM' 
-h 
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A 
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2E„ 



u,. 
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+ 



2E„ 



(4.22a) 
(4.22b) 



(4.23a) 
(4.23b) 



where E^^ are the diagonal matrix elements of the matrix E^ 



E„ 



E^ 



(4.24) 



In this representation, the average pairing gap 
the occupied states, 



is given by the average value of in 



(A) 



(4.25) 



Equations ( |4.23| ) and ( 4.24 ) misleadingly resemble those of the simple BCS theory 
However, in the HFB theory, is not the single-particle energy (i.e., the eigenvalue of h) 
but the diagonal matrix element of h in the canonical basis. Similarly, A^ does not represent 
the pairing gap in the state ip^, and i?^ is not the quasiparticle energy E. However, since 
these quantities define the occupation probabilities, they play a very important role in an 
interpretation of the HFB results, and many intuitive, quantitative, and useful features of 
the BCS theory can be reinterpreted in terms of the canonical representation (cf. Sec. |V C| ). 

In particular, the average values of single-particle p-h and p-p Hamiltonians fulfill the 
following self-consistency equations: 



1 - 



- A ' 

E,. 



A,. 



Eu 



(4.26a) 
(4.26b) 



For a given interaction F^^ and Gfj_u-, Eqs. (|4.26| ) represent a set of nonlinear equations which 
determine and A^. Equations for ( [4.26a| ) and for A^ ( |4.26b|) are coupled by the values 
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of El, ( [4.24| ), which depend on both and A^. However, it is clear that the interaction in 



the p-h channel mainly influences the values of e^, while that in the p-p channel - A 

Unfortunately, Eqs. ( |4.26| ) cannot replace the original HFB equations, because they 
require the knowledge of the canonical basis to determine the F^j, and G^i, matrices. The 
only way to determine the canonical basis is to solve the original HFB equation 



and 

then to diagonalize the density matrix (|4.9a| ). Moreover, solving Eqs. ( [4.26|) ensures that 
only the /i=z/ subset of variational equations ( |4.19a| ) is met, the minimum of energy being 
obtained by solving the whole set (i.e., for all indices n and u). 

The diagonalization of the energy matrix E^^, gives the spectrum of HFB eigenenergies, 

E„: 



HE 



E„U, 



(4.27) 



The matrix U, 
basis 



nfj, represents the unitary transformation from the canonical to the quasiparticle 
Its matrix elements provide the link between the quasiparticle energies En and 
the diagonal matrix elements E^ which define the occupation probabilities, i.e.. 
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n/i I 



(4.28) 



V. COUPLING TO THE POSITIVE-ENERGY STATES 

For weakly bound nuclei one may expect that the particle continuum influences the 
ground-state properties in a significant way. As discussed in Sec. [IV C 2|, the phase space 



corresponding to positive single-particle energies should not be confused with the continuum 
of scattering states which asymptotically behave as plane waves, and are significant for 
genuine scattering phenomena. 



A. Boundary conditions 

Properties of the continuum scattering states are intuitively well understood in terms 
of unpaired single-particle orbits. Shown in Fig. |^ are the self-consistent HF-l-SkP neutron 
single-particle energies in ^^°Sn, e^fj, as functions of the radius i?box of the spherical box in 
which the HF equations are solved. It is assumed that the following boundary condition 
holds for all single-particle wave functions: 

^^(^box) = 0. (5.1) 

For bound single-particle states, e^;'^<0, the effect of increasing i?box beyond 10 fm is in- 
significant. As seen in Fig. ||, the energies of the least bound 3p, 2/, 1/19/2, and 1^13/2 states, 
which form the 82<A^<126 shell, are independent of -Rbox- 

The boundary condition ( |5.1|) leads to a discretization of the continuum by selecting only 
those states which have a node at r=i?box- When i?box increases, the density of the low- 
energy continuum states increases as -Rbox- This effect is very well visible in Fig. Among 
those states whose energies decrease with -Rbox, one may easily distinguish some quasi-bound 
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states, which have energies fairly independent of -Rbox- In Fig. ^ these are the high-£ states 
^11/2, ji3/2, ji5/2, and /ci5/2- However, at some values of i?box they are crossed by, and they 
interact with, the real continuum states (plane waves) of the same quantum numbers, and 
their precise determination is, in practice, very difficult. 

A solution of the HFB equation ( [4. 131) in the spherical box amounts to using the analo- 



gous boundary conditions, 

ME, Rho.) = ME, Rho.) = 0, (5.2) 

for both components of the HFB wave function. As a result, the quasiparticle continuum 
of states with \E\>—X is discretized and becomes more and more dense with increasing 
i?box- However, as discussed in Sec. P^V C[ the density matrix depends only on the localized 
(lower) components of the quasiparticle wave functions and, therefore, is very well stable 
with increasing -Rbox- By the same token, the properties of the canonical-basis states, which 
are the eigenstates of the density matrix, are also asymptotically stable. Of course, the bigger 
the value of -Rbox, the larger is the numerical effort required to solve the HFB equations. 
Consequently, it is important to optimize the value of -Rbox! i-e., to use the smallest box sizes 
which reproduce all interesting physical properties of the system. 

Apart from ours, there are also other possible approaches to solving the HFB eigenprob- 
lem; in particular: (i) the diagonalization in the large harmonic oscillator basis, and (ii) the 
two-step diagonalization. Scheme (i) has been used, e.g., in the HFB+Gogny calculations 
or in the deformed HFB+SkP calculations of Ref. 17911. Its limitations, due to the incorrect 



asymptotics, are discussed in Sec. [V G\ below. In method (ii) one first solves the HF problem 



and then diagonalizes the full HFB Hamiltonian in the HF basis. Such a strategy has been 



suggested in Ref. |H2| and recently adopted in Ref. ||81[]. 



B. Canonical single-particle spectrum 

As discussed in Sec. [IV D|, quantities which determine the p-h properties of the system 



are the canonical energies [Eq. (|4.22a|) ]. The neutron canonical energies in ^^''Sn are 



shown in Fig. [1^ as functions of the box size -Rbox- In this figure, the single-particle index fj, 
is represented by the spherical quantum numbers nij; only the states with occupation prob- 
abilities vl^j>omoi are presented. The canonical states belonging to the shell 82<A^<126 
have negative e„^j's, and they are very close to the HF single-particle energies displayed in 
Fig. ^. They do not depend on the values of -Rbox for -Rbox>10fm. 

At positive values of Cnij, there are several orbitals which do not depend on the box size 
even at -Rbox<15fm. These states correspond to the high-£ quasibound states 211/2, ji3/2, 
ji5/2, and ki5/2, already identified in the HF spectrum of Fig. ^. The values of e„^j for 
these states are only slightly higher than the corresponding values of e^fj. However, these 
quasibound canonical-basis states are not accompanied by the sea of plane-wave scattering 
states (cf. the ji3/2, and the /C15/2 states in Figs. || and |10D. One can thus say that the 
canonical-basis states represent the quasibound states well decoupled from the scattering 
continuum. 

Many other canonical-basis states, especially those with low orbital angular momenta i, 
significantly depend on the box size up to about -Rbox=18fm, and then stabilize. Therefore, 
in all subsequent calculations we use a "safe" value of -Rbox=20fm, unless stated otherwise. 
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Above 20 MeV there appear states with canonical energies fluctuating with -Rbox- These 
states have very small occupation probabilities close to the limiting value of f^£j=0.0001, 
and their determination as eigenstates of the density matrix is prone to large numerical 
uncertainties (see Sec. |I V C 2| ) . One should note that the physical observables are calculated 
directly by using the HFB density matrices, and the above numerical uncertainties do not 
affect the results obtained within the HFB theory. 

As pointed out in Ref. [|7^], the canonical spectrum presented in Fig. can be used 
to analyze the shell effects far from stability. In particular, the size of the A^=126 gap is 
very small (a 2 MeV gap between the 1^13/2 and 4si/2 states), and hence it cannot yield any 
pronounced shell effect (seen, e.g., in the behavior of the two- neutron separation energies, 
Sec. |VI C|) . This shell-gap quenching is not a result of a too small value of the spin-orbit 
splitting. Indeed, a larger spin-orbit strength would push the 1^13/2 level down in energy, 
without affecting the size of the iV=126 shell gap (several negative-parity states are nearby). 
The iV=126 gap, which is equal to about 4 MeV at i?box=10fm, closes up with increasing 
-Rbox due to the several low-£ states whose energies steadily decrease. This effect can be 
attributed to the pairing- induced coupling with the positive-energy states (see Sec. |V G] ). 

In the energy window between and 20 MeV, the density of single-particle canonical 
energies is fairly uniform and no pronounced shell effects are visible. Since the Fermi energy 
must stay at negative values, this region of the phase space cannot be reached. However, 
one may say that the influence of the positive-energy spectrum on the bound states (had 
we analyzed it in terms of, e.g., the Strutinsky averaging) is characterized by a rather 
structureless distribution of states. Above 20 MeV, the occupation probabilities rapidly 
decrease (cf. Table D, and this part of the phase space can safely be disregarded, provided 
one stays in the canonical basis. 



C. Single-quasiparticle spectrum 

The eigenvalues of the HFB equation (^4.13|) (single-quasiparticle energies) carry informa- 
tion on the elementary modes of the system. The lowest single-quasineutron energies E^f^ 
in tin isotopes between A^=50 and A^=126 are shown in Fig. |Tl] (top panel). Apart from 
the magic shell gaps at A^=50 and A^=82, where the single-quasiparticle energies exhibit 
sudden jumps, they depend rather smoothly on neutron number. For a given orbital n^j, 
the minimum of E^f^ is attained in the isotope where the corresponding single-particle state 
is closest to the Fermi energy. Hence, from Fig. [TT| one can infer the order of single-particle 
energies in the beginning of the 50<iV<82 shell as 2^5/2, 3si/2, 2^3/2, 15^7/25 and l/iii/2- 
Similarly, the predicted order at the bottom of the next major shell is 2/7/2, 3^3/2, 3]9i/2, 
2/5/2, 1^9/2, and lii3/2- The order of spherical single-particle states does vary with A^. For 
instance, according to the HFB+SkP calculations of Fig. |ll|, the 1(77/2 shell never becomes 
lowest in energy, as it should have done, had the single-particle energies been A^-independent. 

Noteworthy is the fact that, due to the strong interaction with the low-£ continuum 
(cf. Sec. |V Bp , the 4si/2 excitation becomes lowest at A^>114. Above the 4si/2 state there 
appear several quasiparticle states with excitation energies rapidly decreasing with A^. These 
orbitals represent the low-energy continuum states. They are very close in energy, exhibit 
small spin-orbit splitting, and the lowest of them are the low-£ states: 4pi/2, 4p3/2, 3^3/2, 
and All these features are characteristic of the continuum states Still higher in 
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energy, one may distinguish a similar doublet of the 8/5/2 and 8/7/2 states, as well as the 
2(79/2 state which represents a high-£ resonance. 

The bottom panel of Fig. |Tl] shows similar results for the BCS-like canonical energies 
defined in Eq. ( [4.24|) , and denoted here by -E.^^^. A comparison between -E^™ and E^^^ 
illustrates the fact that the lowest elementary excitations of the nucleus are equally well 
described by both these quantities. Indeed, a general pattern and, in most cases, also the 
values of E^f^ and E^^- are very similar. The differences mainly concern the Si/2 states, 
and also the low-£ states in the continuum, which in the canonical representation appear 
higher in energy (see Table | for the direct comparison for si/2 states). On the other hand, 
the position of the high-£ 2(79/2 resonance is almost identical in both representations. Such 
a similarity supports the supposition (Ref. fT^ and Sec. [VB| ) that the canonical single- 
particle energies, which are the main ingredients of E'^^-, constitute a fair representation of 
single-particle and single-quasiparticle properties of the system. 



D. Relation between canonical and single-quasiparticle wave functions 



As discussed in Sec. [Ill B| , the canonical states constitute a basis in which the 



independent-quasiparticle state has the form of a product of correlated pairs [Eq. ( |8.1ti| )]. 



Therefore, these states can be considered as fundamental building blocks describing the 
pairing correlations in a many-fermion system. On the other hand, the canonical states are 
determined by a solution of the HFB equation - the single-quasiparticle states. 

Since the canonical states constitute an orthonormal ensemble, the lower and upper HFB 
components can be expanded as 

UE^,ra) = Y.A^li,,{ra), (5.8a) 

where 

= / dVE^;(ra)0,(E„,ra) (z = 1,2) (5.4) 

are the associated overlaps. In order to find the relation between A^^^ and A^^^ one can 
employ Eqs. ( p.81|) and ( [4.9|) for the HFB densities. This gives the canonical wave functions 
expressed as linear combinations of the lower HFB components: 



Mr(^)=J2^nlMEn,ra), (5.5a) 

n 

Mr^r) = E4>2(^n,ra). (5.5b) 



One should note that the expansions (|575|) are valid regardless of the fact that the lower 
components (f)2{En,ra) do not constitute an orthogonal ensemble of wave functions. By 
multiplying both sides of Eqs. ( p.5a| ) and ( p.5b|) with iplira) and taking the scalar product, 
one arrives at the orthogonality relations: 
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E^S^it = ^j5.^, (5-6a) 

n 

E^S^i'j* = -WM- (5.6b) 



n 



The above identities express the fact, that both A^^^ and A^^l ^'^^ related to the transfor- 
mation matrix Unu defined in Eq. (|4.27| ) : 

•^n/l ~ l^fJ^nfi ) -^iiii — ~Uijh(niiy (5-7) 



and Eqs. (|5^) reflect the unitarity oiUnv Equations (pTZl ) can be easily derived by inserting 
expansions ( |5.3| ) into the HFB equation ( [4.8| ), and then expressing the matrix E^y ( |4.19b|) 
in its eigensystem ( 4.27| ). 

It is instructive to express the upper HFB component in a form similar to that of Eq. 
(PI): 

ra) = -Y: ^A^^li,,{ra). (5.8) 

For En > —A, the upper component (f)i{En,ra) is the scattering wave function. It can be 
formally expanded in the localized canonical wave functions according to Eq. ( p.8| ), but the 
main contribution comes from the particle-like states with very small values of t;^. Hence, 
this relation is not too useful in practical applications. 



E. Spectral distribution for the canonical-basis wave functions 



In order to discuss the importance of the particle continuum on the structure of canonical 
states, it is interesting to see how a given canonical state is distributed among the single- 
quasiparticle states. For this, it is convenient to rewrite Eq. ( p.5aj ) in the following way: 



0<En<E 

max 



5, 



b2{En,ra). 



(5.9) 



The spectral amplitudes iS„^ define the distribution of the canonical states among the single- 
quasiparticle states. It is important to recall at this point that the sum in Eq. ( ^.9|) represents 
in fact the discrete (£'„<— A) states and the discretized (£'„>— A) continuum states, i.e.. 



E 

0<£;„<-A 



5, 



dn{E)^ME,ra) 

A ^/NE 



(5.10) 



with the spectral amplitudes and iSg^^ pertaining to the discrete and continuous HFB 
spectrum, respectively. 

The spectral amplitudes can be expressed in terms of matrices A^^ or introduced in 
Sec. ITO: 



2 •^"M 



(5.11) 
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We have included in iS„^ the norms Nn of the lower components, Eq. ( [4.11| ). In this way, 



the values of spectral amplitudes measure the influence of quasiparticle states irrespective 
of the overall magnitude of their lower components. 

Before discussing the properties of the spectral amplitudes, let us write down the two 
sum rules: 



^ -^\Snti?j^ -^\Snti\^J^, (5-12) 



i = Ei'5..r^. (5.13) 

The first two sum rules, Eq. (|5.12| ), come from the the unitarity of The last one, Eq. 
( p.l3| ), expresses the condition defining the norm of the lower HFB component. 

In Fig. [l^ are shown the spectral amplitudes for the S1/2 canonical states in ^^°Sn (cf. 



Sees. |1V C 1| and [IV C 21). The phases of the single-quasiparticle wave functions have been 



fixed in such a way that all the amplitudes iS„^ for /i=l are positive (some of these amplitudes 
are too small to be displayed in the figure). This defines the relative phases of the spectral 



amplitudes for /i>l. Then, the positive and negative amplitudes are in Fig. shown by 
bars hashed in opposite directions. Results shown in this figure pertain to the same single- 
quasiparticle and canonical states as those shown in Figs. |^ and 0, respectively, and in Table 

I 

The lowest panel in Fig. |1^ shows that the lsi/2 canonical state (/U=l) is composed mainly 
of two components corresponding to the two deep-hole quasiparticles at £"8=31.64 MeV 
and £5=17.60 MeV. Similarly, the /i=2 and /i=3 canonical states are mixtures of the 
£5=17.60 MeV and £i = 1.54MeV quasiparticles. For all three of these canonical states, 
the diagonal amplitudes dominate. 

Another pattern appears for the positive-energy canonical states; i.e., for /i=4 and 
/i=5. These two canonical states contain large components of the hole-like quasiparticles at 
£5=17.60 MeV and £i=1.54MeV, but in addition, they also acquire large components of 
the particle-type quasiparticles belonging to the continuum. These continuum components 
are centered around 15 and 20 MeV for yU=4 and yU=5, respectively. This illustrates the fact 
that a correct description of the positive-energy canonical states requires solving the HFB 
equation to rather high energies. The widths of the corresponding distributions are rather 
large, which indicates that there is not a single resonance in the particle continuum which 
would alone describe the high-energy si/2 canonical states. This can be well understood by 
recalling that the £=0 resonances have usually very large widths. 

For the drip-line nucleus ^^°Sn, the spectral Si/2 amplitudes are shown in Fig. 13. Sim- 



ilarly to the case of ^^''Sn, the three lowest canonical states for /i=l, 2, and 3 are mainly 
composed of the three hole-like quasiparticles at £9=34.27, £7=22.12, and £3=7.24 MeV 
with dominating diagonal amplitudes. On the other hand, the low-lying positive-energy 
canonical /x=4 state has large and almost equal components coming from the particle-like 
quasiparticles at £i=2.40, £2=4.84, and £4=8.93 MeV. The following yU=5 canonical state 
has dominant amplitudes from the hole-like and particle-like quasiparticles at £3=7.24 and 
8.93 MeV, respectively. One should note that the /i=4 and /i=5 canonical Si/2 states in 
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-^^°Sn have rather large occupation factors as compared to those in ^^"Sn. Both of them 
require including the single-quasiparticle states at least up to 10 MeV. The following /i=6 
state (not shown in the figure) has the occupation probability of fg=0.0003 and the spectral 
amplitudes extending up to 25 MeV. 

The spectral amplitudes for the fi/2 states in ^^''Sn and ^^°Sn are shown in Figs. |1^ and 
|15|, respectively. An interesting situation appears in ^^°Sn where two quasiparticles, one of 
the particle type and another one of the hole type, have rather similar single-quasiparticle 
energies of 17.63 and 18.97 MeV. As a result, the lowest canonical state (/i=l) acquires 
a substantial particle- type quasiparticle component, while both quasiparticles contribute 
almost equally to the /i=3 canonical state. In ^^°Sn, the positive-energy /7/2 canonical 
states (/U=3 and 4) have large amplitudes from the hole-like quasiparticles (contributing 
almost exclusively to the structure of the negative-energy canonical states with /x=l (I/7/2) 
and fi=2 (2/7/2)), as well as from a wide distribution of several particle-type quasiparticles 
extending up to 20 MeV. 

The spectral amplitudes allow also for a determination of the asymptotic properties of 
canonical states. (See Ref. p9| for a discussion of the the asymptotic properties of natural or- 



bits.) The lower components (j)2{En,ra) behave asymptotically as exp{—r y 2m{En — A)//i 
Therefore, as seen from Eq. ( p.lU| ), the asymptotic properties of canonical states 



are governed by the lowest discrete quasiparticle, provided the corresponding spectral am- 
phtude, iSi^, is not equal to zero. However, if such a spectral amplitude is non-zero but very 
small, the corresponding asymptotic behavior will be attained only at very large distances. 
In practice, the lowest discrete quasiparticle dominates the asymptotic behavior only if the 
corresponding spectral amplitude has a significantly large value. For the Si/2 states in ^^°Sn 
(Fig. |12D such a situation occurs for the canonical states with /i=2-5 On the other hand, 
since the value of is very small, the asymptotic behavior of the fi=l canonical state 
is dominated by the hole-like quasiparticle at i?5=17.60 MeV. A similar situation occurs for 
the /7/2 states in i20Sn. Namely, only for the /i=2 canonical state the asymptotic behavior 
is determined by the lowest discrete quasiparticle. 

An entirely different property can occur in drip-line nuclei, where the Fermi energy is 
close to zero and there may be no quasiparticle excitations in the discrete spectrum between 
and —A. In such a situation, shown in Figs. |l^ and |15|, the canonical states are represented 
by superpositions of lower quasiparticle components belonging to the particle continuum. 
Consequently, it is the integral over the lowest continuum quasiparticle states just above 
the E > —A threshold that determines the asymptotic properties of the canonical states. In 
other words, the profile of the level density, dn{E)/dE, around E = —A becomes a crucial 
factor. Good examples of a very strong coupling to the particle continuum are the /i=4 and 
5 canonical Si/2 and /7/2 states in ^^°Sn, where the quasiparticle strength is distributed in 
a very wide energy interval ranging from 1.5 to 20 MeV. On the other hand, the two lowest 
canonical /7/2 states in ^^°Sn can be associated with the two quasiparticle excitations well 
localized in energy (see Fig. 15) and their asymptotics is governed by the energy of the 
lowest quasiparticle. 

An analysis of the spectral distribution, analogous the one presented above, has recently 
been performed for the natural orbits in ^^O determined within the Green's function 
method using the A^A^ interaction. This method accounts for a much more general class of 
correlations as compared to the HFB correlations of the pairing type studied here. However, 
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the general features of the spectral distributions remain essentially the same. Namely, the 
low-occupation-number natural orbits are determined mostly through high-energy contin- 
uum contributions, and large box sizes (15-20 fm) and large single-particle bases (20 states 
per £j-block) have to be used to stabilize the solutions. This is so even if the studied nu- 
cleus (^^O) is /5-stable, well-bound, and light; one can expect that for drip-line nuclei the 
aforementioned features can only be more pronounced. 



F. Asymptotic properties 

In the limit of weak binding, radial dimensions of atomic nuclei increase and it becomes 
exceedingly important to control the radial asymptotics of many-body wave functions, not 
only in reaction studies but also in nuclear structure applications. Figure ^ displays the 
radial dependence of the neutron density p(r) in ^^°Sn calculated with the values of -Rbox 
between 10 and 30 fm. It is seen that, for every value of -Rbox, pij") follows its asymptotic 
behavior up to about -Rbox— 3fm and then falls down to zero as a result of the boundary 
conditions ( ^.21 ). That is, these boundary conditions affect the density only in a narrow 
spherical layer of the thickness equal to about 3fm, while inside this layer p(r) behaves 
independently of the value of -Rbox- Analogous results for the pairing density p(r) are shown 
in Fig. 0. 

At very large distances the asymptotic behavior of the particle density is governed by the 
square of the lower component of the single-quasiparticle wave function corresponding to the 
lowest quasiparticle energy -Emin- Similarly, the asymptotic behavior of the pairing density 
p(r) is determined by the product of the upper and the lower components of quasiparticle 
Ermn- Usiug the asymptotic properties of the HFB wave functions derived in [^,^, one 
obtains: 

, . large r exp(-Xr) I^^A \ 

p{r) — > ; X = 2^2, (5.14a) 

. large r exp(-Xr) . 

p(r) — > ; x = f^i + f^2, (5.14b) 



where 



2m(-E^i„ - A) 2m(E^i„ - A) . 
Ki = y , K2 = y . (5.15) 

In the considered example of ^^'^Sn the calculated values are A=— 1.46MeV and -Einin= 
1.07 MeV (a pi/2 state). Consequently, x^O.TOfm^^ and x~0.49fm~^. In Figs. |16| and |17 



the asymptotic dependencies given by Eq. ( |5.14| ) are shown as shaded lines. One can see that 
for p(r) the asymptotic regime is reached only at distances as large as 25 fm, which means 
that the contributions from other quasiparticle states, and/or from the next-to-leading-order 
terms in the Hankel functions, still influence the particle density at rather large values of r. 
Interestingly, the pairing density approaches the asymptotic limit already at r~10fm. 

A rough estimate of x and x can be obtained by substituting the value of a typical pairing 
gap (A=l MeV) for the lowest quasiparticle energy -Emin- For stable nuclei (A~— 8 MeV) one 
obtains x~1.32fm~^, while for the one-neutron drip nuclei, defined by a vanishing separation 
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energy, S'„^A+A~0, the result is x^0.62fm~^. This difference illustrates the increase in 
the spatial extension of the particle densities when going towards the neutron drip line. On 
the other hand, for the pairing densities the corresponding numbers are x~1.24fm~^ and 
X=x/2~0.31 fm~^. Therefore, in stable nuclei both types of densities have rather similar 
asymptotic behavior, while in drip-line nuclei the pairing densities have much longer tails. 

In this context, it is instructive to recall the discussion from Sec. [Ill C| regarding the 
probabilistic interpretation of the HFB densities. The probability Vi{x) {V2{x)) of finding 
a particle or a pair of particles at r=x is proportional to p(x) or p'^{x) + p^{x), respectively. 
Consequently, in stable nuclei V2{x) decays much faster than Vi{x) at large distances. This 
is not true for drip-line nuclei, where the asymptotics of Vi{x) and V2{x) is the same. 

As discussed above, static pairing correlations can influence dramatically the asymptotic 
behavior of density distributions in drip-line nuclei. In addition, a significant modification of 
the density tails comes from the dynamical coupling to collective modes through the particle 
continuum. Such a coupling can be treated in terms of the continuum QRPA and has been 
shown to be very important for light systems |8^j8^ . An analysis of the asymptotic behavior 
of the particle density p(r) has recently been performed ||8^ by finding the exact solutions 
for weakly bound two particles interacting through a contact force. In that study, the role 
of one-particle resonant states on the density asymptotics has been discussed. 



G. Pairing coupling to positive-energy states 



As illustrated in Sec. |V A| , the density of the scattering continuum states increases with 
i?box- In the limit of very large values of Rhox, the set of discretized continuum states 
can be considered as a fair approximation of the real continuum, and the sums over the 
positive-energy states can correctly represent integrals over the continuous energy variable. 
Therefore, we may consider this limit in order to study the dynamical coupling between 
the bound single-particle states and the positive-energy states. In the language of pairing 
correlations, one may think of this coupling in terms of a virtual scattering of pairs of 
fermions from the bound states to positive-energy states, and back. Such a pair scattering 
gives rise to the additional pairing energy to the ground-state energy. 



To illustrate the stability of results with increasing box size, in Fig. |I8| we show the 
neutron p-p potentials U{r) in ^^°Sn and ^^^Sn calculated in the HFB-|-SkP model for several 
values of -Rbox- In these two nuclei, the values of f/(r) do not change when Rhox is larger 
than 20 and 22 fm, respectively, but at smaller values of -Rbox, one observes significant 
variations. A rather unexpected result of this analysis is that the overall magnitude of 
pairing correlations, represented by the average pairing gap (A), decreases with increasing 
i?box- This occurs in spite of the fact that the actual density of scattering states dramatically 
increases with increasing i?box- 

This effect can be understood by noting that the pairing correlations produced by a 
density-dependent p-p interaction (and hence for the SkP force used here) are concentrated 
at the nuclear surface; i.e., at a fixed location in space. For small values of -Rbox, the 
boundary conditions ( |5.2|) have a tendency to push the continuum wave functions towards 
smaller distances, and into the surface region. This increases the magnitude of pairing 
correlations. On the other hand, with increasing -Rbox, the scattering states spread out 
uniformly outside the nucleus and effectively leave the surface region. Hence (A) decreases. 
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As a consequence, with increasing i?box the self-consistent attractive pairing potential f/(r) 
decreases in magnitude and significantly spreads out towards large distances. 

The importance of allowing the pairing interaction to couple properly to the particle 
continuum is illustrated in Fig. 19, where the neutron rms radius, the average pairing gap, 
and the Fermi energy are shown as functions of -Rbox- The two upper plots confirm that a 
stability of results is attained beyond 20 or 22 fm, while the bottom plot indicates that the 
pairing coupling to the positive-energy states can be a decisive factor infiuencing the nuclear 
binding. Indeed, below i?box— 20fm the nucleus ^'^^Sn is unbound, and it becomes bound only 
when its ground state is allowed to gain an additional binding from the pairing correlations 
at large distances. This indicates that, for the surface-type pairing interaction, one has to 
consider a rather dense particle continuum before the pairing coupling to positive-energy 
states is exhausted. (For a similar discussion in a schematic model see Ref. [Q. There, it 
has been pointed out that because of strong coupling to the continuum, A is significantly 
lowered in the case of surface pairing as compared to the case of volume pairing.) 

Since, for the Gogny interaction, the HFB equations are solved by expansion in the 
harmonic oscillator basis, one can test the coupling to the positive-energy states by in- 
creasing the number A^^sh of the oscillator shells used in the basis. In practice, calculations 
must be restricted to Nsh^^O, which allows one to describe the wave functions up to about 



Rj^^^c:i^j2Nshh/mujQ, where is the frequency of the harmonic oscillator For Nsh=20 
this corresponds to about i?max=14fm. 

Figure compares the asymptotic behavior of the neutron particle densities in three 
neutron-rich tin isotopes calculated in the spatial coordinates (SkP) or in the harmonic- 
oscillator basis (DIS). In the former case one obtains a clean region of the asymptotic de- 
pendence governed by Eq. (|5.14aj) , which around r=18 fm is perturbed by the box boundary 
conditions (|5.2|) at i?box=20 fm. In the latter case, the region of proper asymptotic behav- 
ior becomes perturbed by the exp{—mujQr'^/h) dependence characteristic of the harmonic- 
oscillator-basis wave functions. The ujq values, obtained by minimizing the total energy for 
the iVsh=17 basis, are equal to 13.4, 6.6, and 6.3 MeV and '^'^^Sn, respectively. 

Due to this, a study of the continuum infiuence using such a basis can be performed only up 
to densities of scattering states corresponding to about -Rbox=14fm in the heavier isotopes 
and only i?box=10fm in ^^^Sn, as can be seen in Fig. ^ Let us note, however, that the 
neutron densities beyond r=10fm are typically smaller than 10~^fm~^, which explains the 
stability of the HFB calculations with increasing size of the basis. 

This is illustrated in Fig. ^ which is analogous to the similar study presented for the 
SkP interaction in Fig. |T^. Here, for each value of Ash and for each nucleus, the value 
of Uq was optimized so as to minimize the total energy. As can be seen, one obtains a 
nice stability of results by using A^sh=17. This test corresponds to testing the coordinate- 
representation solutions (Fig. |19]) in the range of box sizes between 12 fm<i?box<14fm. In 
this rather narrow region, the SkP results are not stable because of the dominant surface- 
type character of its pairing interaction. Since the p-p Gogny interaction is more of the 
volume type (Sec. |111 C Ij) it requires much smaller distances to saturate. 
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H. BCS approximation 



When inspecting Fig. it is obvious that by applying the BCS approximation to the 
state-independent pairing force and by allowing the BCS-type pairing correlations to develop 
in such a dense spectrum, the result can be disastrous. The seniority force gives rise to the 



non-localized pairing field fS^ , 



hBCs{ra,r'a') = -ABcs'5(r- - r')6^.', (5.16) 

i.e., to a constant pairing gap, identical for all states. The high density of single-particle 
states in the particle continuum immediately results in an unrealistic increase of BCS pairing 
correlations []T1|. One may, in principle, artificially readjust the pairing strength constant to 
avoid such an increase, but then the predictive power of the approach is lost and, moreover, 
the spatial asymptotic properties of the solutions are still going to be incorrect. 

To illustrate the latter point. Fig. ^ (top panel) shows the neutron densities in ^^°Sn 
calculated for several values of -Rbox within the HF-I-BCS approximation. In order to avoid 
the increase of pairing correlations with increasing density of states, the calculations have 
been performed by fixing the values of the pairing gap. For every box size -Rbox? the value of 
Abcs has been set equal to the average pairing gap (A) obtained within the HFB method. 
The corresponding (A) values are quoted in Fig. 



It is not too surprising to see that the asymptotic behavior of the density calculated 
in the HF+BCS+(A) method (top panel) is entirely different than that shown in Fig. 
Due to a nonzero occupation probability of quasibound states, there appears an unphysical 
gas of neutrons surrounding the nucleus. In Fig. ^ this gas has a constant density of 
p~6x 10~^ fm~'^, independent of -Rbox- This result means that an external pressure would 
have been necessary to keep the neutrons inside the box. Namely, had the box boundary 
condition been released, one would have observed a stream of neutrons escaping the nucleus. 
This is a completely artificial (and unwanted) feature of the BCS approximation, because 
for a negative value of the Fermi energy, neutrons cannot be emitted. 

In the above example the density of the neutron gas at i?box=25 fm corresponds to about 
4 neutrons uniformly distributed in the sphere of i?=i?box- Needless to say, by increasing the 
box radius, the number neutrons in the gas grows at the expense of the number of neutrons 
constituting the nucleus in the center of the box. Since the total average number of neutrons 
is conserved, by changing i?box one actually performs an unphysical study of different nuclei, 
surrounded by a neutron gas of a fixed density. Another consequence of the presence of a gas 
of particles is that the rms nuclear radius cannot be calculated in the BCS theory, because 
the results strongly depend on the box size (see discussion in Refs. [p^ , pJj] ). 



It has been suggested in the literature [gg] that the above deficiencies of the BCS approx- 
imation can be cured by applying to them the state-dependent-pairing-gap version, where 
the pairing gap is calculated for every single-particle state using an interaction which is not 
of the seniority type. (The corresponding BCS equations resemble the canonical-basis rela- 
tions (|4.26| ).) In such an approach one hopes that the majority of continuum states would 
neither contribute to the pairing field (e.g., because of their very different spatial charac- 
ter) nor result in the appearance of the unphysical gas. This conjecture is tested in Fig. ^ 
(middle and bottom panel) where the neutron densities obtained within the state-dependent 
version of the BCS approximation using the SkP'^ and the SkP interactions are presented. 
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It is seen that a reduced coupling of some continuum states to the pairing field does indeed 
decrease the gas density, however, the asymptotic behavior of the density is still incorrect. 

In the above plots, the shaded lines represent the asymptotic behavior given by Eq. 
( ^.14a| ) assuming E^i^=0, i.e., that of a single-particle state at the Fermi energy. It is seen 
that a surplus density above this asymptotic limit appears at large distances. However, 
the deficiencies of the state-dependent BCS approximation, as used for example in Refs. 
5B| , pB|J57| , are certainly less acute than those of the seniority-pairing BCS. For example, in 



this type of approach one may probably calculate radii of nuclei much nearer to the drip 
line. 

It is clear that the neutron gas appears in the BCS solutions because of the nonzero 
occupation probabilities of scattering states. Therefore, one may think that excluding the 
scattering states from the pairing phase space could be a decisive solution to the problem. 
However, for drip-line nuclei, where the Fermi energy is by definition close to zero, the re- 
maining phase space would then be small, and this would lead to an artificial quenching 
of pairing correlations. Moreover, even if the density obtained in such method would van- 
ish asymptotically, the corresponding factor x would not be governed by A— A~2MeV, as 
discussed in Sec. |VF| , but by the single-particle energy, e~0, of the highest-energy single- 
particle state considered in BCS calculations. This again would lead to densities vanishing 
at much slower pace than it is required by the HFB theory. 

VI. PHYSICAL OBSERVABLES FAR FROM STABILITY 

In this section discussed are some experimental consequences of the HFB theory, partic- 
ularly important for weakly bound nuclei. 

A. Pairing gaps 

Pairing gaps are p-p analogs of single-particle energies. They carry the information about 
the energies of non-collective excitations, level occupations, odd-even mass differences, and 
other observables. The average neutron canonical pairing gaps ( [4.22b| ) are shown in Figs. 
^ (^^°Sn) and ^ (^^°Sn) as functions of the canonical single-particle energies ( [4.22aj) . 



As seen in the middle part of Fig. ^ pairing gaps obtained with the volume-type pairing 
interaction exhibit very weak configuration dependence. In ^^°Sn they decrease slightly 
with but remain confined between 1.0 and 1.5 MeV. In general, the values of for the 
si/2 states are slightly larger than for other orbitals, which is again related to the volume 
character of volume delta interaction. 



The results presented in the bottom part of Fig. ^ nicely illustrate the surface character 
of the SkP pairing interaction. Indeed, here the pairing gaps increase from 0.5 MeV (deep- 
hole states) to about 1.25-1.5 MeV when the single-particle energies increase towards the 
Fermi energy, and then they decrease again to about 1.0 MeV for positive single-particle 
energies. This is related to the fact that orbitals near the Fermi level are concentrated in 
the surface region. 

Still another type of behavior is obtained for the finite range Gogny interaction (top 



part of Fig. ^S]). Here, the pairing gaps decrease steadily with single-particle energy. In 
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-^^°Sn the values of decrease from about 2.5 MeV for deep-hole states to about 0.75 MeV 
for positive-energy states. (A similar energy dependence of pairing gaps was obtained in 
the BCS calculations of Ref. [|1^ with the renormalized Paris potential.) Interestingly, the 
values obtained for the high-£, j=(i—\ orbitals (antiparallel L—S coupling) are significantly 
larger than those for other orbitals. The different ranges of values for SkP and DIS in 
Fig. ^ reflects the different effective masses in both models. A rather low effective mass 
in DIS, m*/m=0.70, gives rise a reduced level density and a more bound lsi/2 ground 
state as compared with the SkP model (m*/m=l). In fact, due to the non-local exchange 
contributions to the p-h mean field (Appendix ^), the lsi/2 state in the Gogny model has 
the canonical energy lower than the bottom of the local potential well, shown in Fig. ^. 

In ^^°Sn, the HFB+DIS pairing gaps at the Fermi energy are of the order of 1.75 MeV, 
which slightly overestimates the values corresponding to the odd-even mass staggering in 
this region. However, one should bear in mind that the pairing gaps at the Fermi energy are 
rather rough approximations to the odd-even mass difference. A more accurate description 
can be obtained by performing blocked HFB calculations for odd-mass isotopes. In the 
vicinity of ^^°Sn this method yields the odd-even mass staggering of 1.6 MeV ||5^ for the DIS 
interaction and of 1.3 MeV [0] for the SkP interaction. Another contribution to the odd- 
even mass difference comes from the coupling to the low-lying collective modes. Therefore, 
the DIS parameters have been adjusted to give the pairing gap in tin to be 0.3 MeV 
larger than the experimental one. On the other hand, such a margin has not been taken 
into account for the SkP and SkP*^ forces. Clearly, a detailed comparison of the values of 
pairing gaps for the interactions discussed in Fig. ^is delicate. Much more information can 
actually be derived from the comparison of their dependence on the single-particle energies, 
which is markedly different. 

The general pattern of remains very similar when going to the neutron-rich nucleus 
^^°Sn (Fig. 1^). In particular, the magnitude of the average pairing gap in deep- hole states 
depends strongly on the range and density dependence of pairing interaction. 

Figures ^ shows the average neutron pairing gaps (Eqs. ( [4. 61) and ( |4.25| )) for SkP, SIII'', 
and DIS interactions. The large values of (A) obtained in HFB+DIS can be explained by: 
(i) an overall larger magnitude of pairing correlations in tin nuclei, and (ii) strong pairing 
correlations in deep-hole states which strongly contribute to the average, Eq. (|4.25| ). It is 
to be noted, however, that despite stronger pairing in DIS, the HFB+DIS pairing gaps 
vanish at A^=126 (near the two-neutron drip line), in contrast to the HFB+SkP result. This 
difference may be traced back to a much larger continuum phase space taken into account in 
our HFB+SkP calculations (Sec. [VG|) which are performed in the coordinate representation, 
and to a larger A^=126 shell gap (4.2 MeV in ^^^Sn) obtained with DIS. (The increase of 
proton pairing gaps when approaching the proton drip line has been calculated previously in 
Ref. with the HFB+SkP model and explained in a similar way.) The disappearance of 
the neutron pairing at A^=126 in the HFB+SIII'' model is partly due to the volume character 
of h (a weaker coupling to the particle continuum) and partly due to a larger iV=126 shell 
gap 
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B. Shell effects 



As discussed in Sec. [IV A| , diffused nucleonic densities and very strong, surface-peak;ed, 
pairing fields obtained with the density-dependent pairing interaction are expected to lead 
to very shallow single-particle potentials in drip-line nuclei. Because of a very diffuse surface 
(no flat bottom), the resulting single-particle spectrum resembles that of a harmonic oscil- 
lator with a spin-orbit term (but with a weakened £^ term) |7^. Schematically, this effect 
is illustrated in the left panel of Fig. E6l By comparing with the situation characteristic of 



stable nuclei (right panel of Fig. ^6]), new shell structure emerges with a more uniform dis- 
tribution of normal-parity orbits, and the unique-parity intruder orbit which reverts towards 
its parent shell. Such a new shell structure, with no pronounced shell gaps, would give rise 



to different kinds of collective phenomena [ff^ 



The effect of the weakening of shell effects in drip-line nuclei, first mentioned in the 
astrophysical context [^, was further investigated in Refs. [^,[7^,Q. First analyses of its 
consequences for the nucleosynthesis have also been performed Microscopically, it 

can be explained by: (i) the changes in the mean field itself due to weak binding (see above), 
and (ii) a strong pairing-induced coupling between bound orbitals and the low-£ continuum. 

C. Separation energies 

Weakening of shell effects with neutron number manifests itself in the behavior of two- 
neutron separation energies. This is illustrated in Fig. ^ which displays the two-neutron 
separation energies for the A^=80, 82, 84, and 86 spherical even-even isotones. The large 
N=82 magic gap, clearly seen in the nuclei close to the stability valley and to the proton- 
drip line, gradually closes down when approaching the neutron drip line. A similar effect 
is seen in the [Z, N) map of the spherical two-neutron separation energies for the particle- 



bound even-even nuclei calculated in the HFB-(-SkP model (Fig. |2^). Namely, the neutron 
magic gaps A^=20, 28, 50, 82, and 126, clearly seen as cliffs in the 5*2™ surface, disappear for 
neutron-rich systems. 

The gradual disappearance of the neutron shell structure with N is not a generic property 



of all effective interactions. As seen in the plot of S2n and Xn for the tin isotopes (Fig. ^9]) 
this effect is seen in the SkP and SkP'^'' models, and, to some degree, also in the SkP*^ model. 
(A weak irregularity at A^=126 reflects the weaker coupling to continuum for the volume 
pairing |3^.) The strong shell effect seen in the SIII and SkM* results has been discussed 
in Ref. [Q; it can be attributed to the low effective mass in these forces. The result of the 
DIS model, both for and X^, is close to that of the SkP^ model. It is interesting to 
point out that the QLM calculations of Ref. (with m*/m = 1) for the Sn isotopes yield 
very similar results to those of HFB-|-SkP. 

The very neutron-rich nuclei, as those shown in Fig. cannot be reached experimentally 
under present laboratory conditions. On the other hand, these systems are the building 
blocks of the astrophysical r-process; their separation energies, decay rates, and cross sections 
are the basic quantities determining the results of nuclear reaction network calculations. 
Consequently, one can learn about properties of very neutron-rich systems by studying 
element abundances The recent r-process network calculations [|13|, based on several 
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mass formulae, indicate a quenching of the shell effect at N=82 in accordance with the 
results of HFB+SkP model. 



D. Deep hole states 



Pairing interaction between bound orbitals and particle continuum is partly responsible 
for the appearance of particle widths of deep-hole states and the term-repulsion phenomenon 
(strong repulsion between single-particle levels) [p3| , p4| . In the DWBA and for the local 
pairing field U the particle width is given by 



2n 



(6.1) 



Here, (fi{r) is the HF wave function of the bound deep- hole state i with the single-particle 
energy X — Ei in the absence of pairing, while ife{f ) is the HF wave function of the unbound 
state with the energy \ + Ei. 

Equation (|6.1|) is obtained by assuming that the p-p field of the HFB Hamiltonian can 
be treated perturbatively. A more consistent way would be to estimate Fj based on self- 
consistent HFB solutions containing pairing correlations. The proper formulation of the 
nonperturbative HFB-based theory of deep hole states and one-particle transfer process still 
needs to be developed. 

As discussed in Ref. [34], Fj is sensitive to the type of the pairing force. In general, the 
widths are larger for surface pairing than for volume pairing. However, the result for an 
individual state strongly depends on its angular momentum and excitation energy. 

Experimentally, total widths of deep hole states, Ftot, are of the order of MeV's, (see, 
e.g., Refs. |9^-p5|). That is, the partial width ( |6.1|) , of the order of 10-100 keV, constitutes 
an extremely small fraction of Ftot- Consequently, the experimental determination of Fj 
alone is very unlikely. 



E. Pair transfer form factors 

There are many interesting aspects of physics of unstable nuclei which are related to 
reaction mechanism studies: weak binding, large spatial dimensions, skins (see, e.g., Refs. 



J96|j97|| ). Below, we discuss some consequences of surface-peaked pairing fields for pair 
transfer studies. 

An experimental observable that may probe the character of the pairing field is the 
pair transfer form factor, directly related to the pairing density p. The difference in the 
asymptotic behavior of single-particle density p and pair density p in a weakly bound system 



see Sees. [Ill (J 1| and |V b'\) can be probed by comparing the energy dependence of one-particle 



and pair-transfer cross sections. Such measurements, when performed for both stable and 
neutron-rich nuclei, can shed some light on the asymptotic properties of HFB densities; 
hence on the character of pairing field. 

Figure RO displays the pair transfer form factors r^p(r) calculated in ^^*^Sn, ^^°Sn, and 
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Sn with the SkP interaction. These microscopic results are compared with the macroscopic 
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form factors r'^Sp{r) which are determined by using the derivative of the particle density 
with respect to the neutron number: 

where -Epair is given by Eq. ( |4.5| ). This expression can be motivated by the fact that only 
the orbitals near the Fermi surface make significant contributions to the pair density. In the 



BCS theory, the normalization constant in Sp{r) is usually chosen as A/G=— impair/ A. 
Here, we use neither the BCS approximation nor the constant pairing strength G. Therefore, 
the normalization — -Epair/ (A) is employed. The derivative in Eq. ( |6.2| ) is calculated from 
the finite difference between the self- consistent results for the HFB vacuum corresponding 
to particle numbers and A^— 1. In these calculations, in order to explore the smooth 

dependence on the particle number N, the odd-average-particle-number vacua have been 
calculated without using the blocking approximation. It should be mentioned at this point 



that the further approximation [p^ , p.OCI|| of the derivative dp{r)/dN by the spatial derivative 



dp(r)/dr is not justified, because the volume-conservation condition is not valid for the 
neutron density distribution (see Fig. |I|). 

The pair transfer form factors in Fig. ^ clearly show that this process has a predom- 
inantly surface character. The macroscopic form factors have smaller widths and higher 
maxima than the microscopic ones. On the other hand, they are smaller in the interior of 
the nucleus as well as in the asymptotic region. In /3-stable nuclei the macroscopic approxi- 
mation works fairly well, while in the drip-line nuclei the differences between the two form 
factors are markedly larger. In general, the corresponding differences are much larger than 
those obtained within the BCS and the particle-number-projected BCS approaches for the 
seniority interaction [|101|| . 

A comparison of the results obtained for different isotopes conspicuously shows a signif- 
icant increase in the pair transfer form factors in the outer regions of drip-line nuclei. In 
^^''Sn, the form factors vanish around 9fm, while in ^^°Sn and ^^^Sn they extend to much 
larger distances. This effect is particularly pronounced for the microscopic pair transfer form 
factors. 



F. Other observables 



The importance of the HFB treatment for calculations of nuclear radii has been discussed 
in several papers ||5^ , [79|j3T| , p!l|] . As mentioned in Sec. y, odd-even staggering of rms charge 
radii is one of the best experimental indicators of the density-dependent pairing. The proper 
treatment of the pairing effect on radii is especially important for weakly bound systems 
which exhibit halo or skin effects [^,^,|ll|] (cf. discussion in Sec. |V G|) . 

Apart from the information on the nuclear rms radii, one may also gain some experimen- 
tal insight into the ratios of neutron and proton densities at large distances from the center 
of nucleus [ [102|jl03|| . This is possible due to experiments on antiproton annihilation from 
atomic orbits, which leads to different reaction products depending on whether the process 
involves a proton or a neutron. 

The role of deformation in neutron drip-line nuclei still needs to be investigated. One 
can anticipate that due to: (i) very diffused surfaces, and (ii) strong pairing correlations, the 
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geometric concept of collective deformation (defined as a deviation of miclear surface from 
sphericity) should be revisited. In this context, the symmetry-unrestricted HFB calculations 
in coordinate space are called for. 

VII. SUMMARY AND CONCLUSIONS 

The advent of radioactive nuclear beams provides many exciting opportunities to create 
and study unstable nuclei far from the j3 stability valley. One of the unexplored areas 
far from stability is physics of nuclear pairing in weakly bound nuclei, especially near the 
neutron drip line. Contrary to the situation characteristic of stable nuclei, the coupling 
between the p-h field and the p-p field in nuclei with extreme N/Z ratios is dramatic; i.e., 
no longer can pairing be treated as a residual interaction. 

The main objective of this study was to perform a detailed analysis of various facets of 
pairing fields in atomic nuclei. The first part contains the comprehensive summary of the 
HFB formalism, with particular attention on the physical interpretation of the underlying 
densities and fields. Very little is known about the p-p component of the nuclear effective 
interaction; its structure is of considerable importance not only for nuclear physics but also 
for nuclear astrophysics and cosmology. Therefore, the second part of this work focuses 
on the differences between various pairing interactions. In particular, the role of density 
dependence and finite range of the p-p force has been illuminated, and the importance of 
the coupling to the particle continuum has been emphasized. Finally, the third part of 
our study relates the theoretical formalism to experimental observables; i.e., energy spectra, 
masses, radii, and pair transfer form factors. It is demonstrated that these observables carry 
invaluable information that can pin down many basic questions regarding the effective NN 
force, and its pairing component in particular. It should be stressed, however, that in order 
to see clearly some of the predicted effects, the excursion far from the valley of /3-stability 
is necessary. 

The analysis presented in this paper should be viewed as a useful starting point for future 
investigations. One of them is the coupling between collective surface modes (e.g., deforma- 
tion) and pairing fields in weakly bound nuclei. Another interesting avenue of explorations 
is the role of dynamics; e.g., the importance of the particle number conservation and the 
coupling to pair vibrations. A fascinating and difficult research program is the microscopic 
description of excited states, especially those lying above the particle emission threshold, for 
which the boundary conditions used in this study (an impenetrable box) have to be modified 
to account explicitly for outgoing waves. We are only beginning to explore many unusual 
aspects of the nuclear many-body problem offered by systems with extreme N/Z ratios. 
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APPENDIX A: THE P-H AND P-P MEAN-FIELD HAMILTONIANS FOR A 
LOCAL TWO-BODY FINITE-RANGE GOGNY INTERACTION 



The Gogny force 0,0] is composed of the central, spin-orbit, density-dependent, and 
Coulomb interactions. The spin-orbit and density-dependent terms have zero-range, and 
their contributions to the p-h and p-p mean fields are identical to those of the Skyrme inter- 



action. The corresponding expressions can be found in several papers; e.g., Refs. ||104| , |52 



and will not be repeated here. In the following we only consider the central finite-range and 
Coulomb terms. The central components read 



2 (-r-r'. 



Ken = 2^e iW, + B,P^-H,P,-M^P^P,), (Al) 

where P^ and Pr are the exchange operators for spin and isospin variables, respectively. 
This interaction is local; i.e., it should be multiplied by 6{ri — r[)6{r2 — r'2) before it is 
inserted in the integrals ( |4.3| ) defining the mean fields. Moreover, it should also be multiplied 
by the antisymmetrizing operator (1 — PrPaPr), where Pr is the exchange operator for 
space variables. One usually calls the term involving Pr the exchange term, while the term 
involving no space exchange is called the direct term. 

The space, spin, and isospin variables are denoted by r, cr=±|, and r=±|, respectively. 
The parameters fij, Wj, Bj, Hj, and Mj, belong to the set called DIS ||105|| which has been 
used in this paper. Since the expressions given by the j=l and 2 components are identical, 
in what follows we drop the index j to increase the legibility of the formulae. 



1. Contribution of the central direct interaction to the p-h mean field 

Since the interaction (Al) is local, the direct term gives the p-h mean field ( |4.3a| ) which 
is also local, i.e.. 



+ 6{r - r')6^^^, f d'r.e ^(fi - M^.^Jp^H^a, n - a), 

where p'^(r) is the density of nucleons (|3.30a|) of type r. 

Assuming that we consider only the states which are even with respect to the time 
reversal, the density matrix ( |3.17a|) obeys the relation ( p.l8a| ). Consequently, the densities 
p'^{ra, rcr) for cr=±| are equal to |p'^(r), and the densities p^ira, r — a) vanish. Therefore, 



(A2) 
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the term in ([X^), which is proportional to S^-a'-, vanishes, and the contribution of the direct 
term to the p-h mean field is the local, spin-independent potential: 



(A3) 



where 



U{r) 



(A4) 



One should note that due to the locality of the interaction, the direct term depends only on 
the local densities. 



2. Contribution of the central exchange interaction to the p-h mean field 



Due to the locality of the interaction, the contribution of the exchange term to the p-h 
mean field involves no integration: 



+S^^^, {H - W5rr, )p^' (ra, r' - a) 



(A5) 



Here the time-reversal symmetry does not bring any simplification. However, a simpler 
formula is obtained in cases where p^^ {ra, r'a') is real. It follows from Eq. ( |3.18a| ) that the 
densities p'^^(rcr, r'a) are equal to | J2a P^^{f<^, f'<^), which finally leads to 



r:,,(ra, r'a) = E + H/2 - {B + W/2)6rrj] E P"' (^^i' ^'^0' (A6) 

3"! 

(A7) 



Tl^^ira, r'-a)=e-^J2iH- 1^5.n)p" (ra, r' - a). 



3. Contribution of the Coulomb interaction to the p-h mean field 

Derivation of the direct and exchange Coulomb fields is similar to the one of the finite 
range term (^) with several additional simplifications. When the nuclear state is time- 



reversal invariant, one obtains the following contributions to the proton p-h mean field in 
terms of the proton densities: 



2 

rLui-dir(^^, r'^') = Sir - r')5^„. j dVi^-^//(rO, (A8) 
rLui-exc(^^, r'a') = rr^/(ra, r'a'). (A9) 
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4. Contribution of the central interaction to the p-p mean field 



The general form of the pairing field is given by Eq. ( [4.3 b| ). In this case the direct and the 
exchange contributions are equal. For the local central force the total contributions 

to the p-p mean field have the form: 

h^{ra,r'a) = e~^^^l^[{W - H)~p'{ra,r'a) - {B - M)~p' {r' a,ra)l (AlO) 
h7{ra, r' - a) = e~^^^^ {W + B-H - M)p^{ra, r' - a). (All) 

Again it is to be noted that due to the locality of the interaction, the corresponding p-p mean 
fields do not involve any integration but are proportional to the pairing density matrices. 
In the case considered in this study (time-even densities), the contribution (|A11|) vanishes. 

Since the exchange parameter of the zero-range density-dependent term of the Gogny DIS 
interaction is fixed at Xq=1, this term does not contribute to the p-p mean field. Moreover, 
the spin-orbit and Coulomb terms usually give small contributions as compared to those of 
the central force (|AT|). 



5. Numerical methods used for the calculation of the mean fields 

Computation of the exchange p-h mean fields, Eqs. (|AB|), (|A^), and (|A^), and the pairing 
fields, Eqs. (|A10|) and ( |A11| ), is straightforward. It only requires the knowledge of the spatial 
spin-dependent non-local particle, p'^{ra,r'a'), and pairing {ra,r' a') densities. 

Computation of the direct p-h mean field, Eq. ([A3|) , and the direct Coulomb mean field, 
Eq. ( |A8|) , is more complicated since it requires the evaluation of three-dimensional integrals 
of the form: 

I,{r) = I d'r'e-^ p(r'), (A12) 
Ic{r) = f dV^- p(r'). (A13) 

In order to compute Ifj,{r) of Eq. ([A12| ) we note that for the single-particle wave functions 
expanded in the harmonic oscillator basis, the local density p(r) is the product of a Gaussian 
factor and of a polynomial in the spatial coordinates Xi, X2, and 



p{r) = exp 



Xk\^ 



k=l ^^kJ 



P(X1,X2,X3), (A14) 



where 6i, 62 and 63 are the HO lengths of the basis states. Consequently, this integral can 
be evaluated exactly using the Gauss-Hermite quadrature. 

The computation of the Coulomb integral ( |A13| ) is more difficult due to the infinite range 
of the Coulomb force. The method we have used consists of expressing the Coulomb force 
as a sum of Gaussians: 

1 2 /■°° d/i .iii-V 

1 7\ = ^ / (A15 

\r — r'\ ^Tx Jo p 
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Then one obtains 



Hr) = %I,{r) (A16) 

y/H Jo /i^ 



where Ifj,{r) is given by Eq. (|A12| ). In order to perform the remaining one-dimensional 



integration, the variable is changed to 



e = &/\/&'+/i^ (A17) 

where b is the largest of the three harmonic-oscillator lengths hi, &2 and 63. This change of 
variable is very convenient, since then the range of integration becomes [0, 1]. The integral 
in Eq. ( |A16| ) can be very accurately computed using the Gauss-Legendre quadrature. 



APPENDIX B: THE ENERGY CUT-OFF 

Calculations which are based on the schematic pairing interaction, or on the contact force 
[Eq. (|2.2| )] require a finite space of states in the p-p channel. For such interactions, when this 
space is increased, the pairing energy diverges for any fixed strength of the interaction. This 
divergence is a well-known effect related to the fact that for the contact interactions the 
matrix elements do not (or too slowly) decrease with the excitation energy. This is not a 
case for finite-range interactions, such as Gogny, for which the pairing energy converges to 
a finite value. 

Since it is considerably easier to use the zero-range interactions than the finite-range 
interactions, one applies the former ones in a limited configuration space determined by 
a cut-off in the single-particle energy or in the single-quasiparticle energy. This can be 
understood as a phenomenological introduction of the finite range |^ . There are two other 
arguments in favor of such a procedure. First, the scattering of particles in the nuclear 
medium at very high energies (or at very small distances) is very little known, and the 
particular form offered by any phenomenological finite-range force is very uncertain. Second, 
the single-particle wave functions are primarily determined by the p-h channel of interaction, 
and they, in general, spread throughout distances which are much larger than the range of 
the p-p interaction. Therefore, physical differences between the zero- and short-range p-p 
forces cannot be expected to be very pronounced. 

Within the BCS approximation, and assuming a constant density of the single-particle 



states at large energies, one can derive [p3| , |106|] a prescription to renormalize the strength 



of the p-p interaction in such a way that the pairing gap A does not depend on the energy 
cut-off. Suppose that the single-particle states with energies —u<e—X<eu are used to solve 
the BCS equations for the force of strength Vq. Then, within the specified approximations, 
the following relation holds: 

where Cq— 300 MeVfm^ is a constant inversely proportional to the density of single-particle 
states near the Fermi energy. In other words, for given values of Cq and A, Eq. ( pi| ) gives 
values of Vq for any other choice of the cut-off energies ej and e'^. 
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Since in the present study we use the HFB method instead of the BCS approximation, 
and since the density of states can hardly be considered to be independent of energy (actually 
for fixed i?box it increases as a/c), the formula ( [BTD cannot be directly used. However, the 
question as to what extent the pairing strengths can be renormalized for a zero-range p-p 
interaction can be addressed by analyzing the numerical solutions of the HFB equations. 

Figure ^ (top panel) shows the neutron pairing energies -Epair [Eq- (Eil)] calculated for 
the SkP'' interaction which uses the contact p-p interaction ( ^.2[ ) with Vo=— 160 MeVfm^. 
It should be recalled at this point that for all coordinate-space HFB calculations presented 
in this study, the cut-off energy, -Emax, depends on the quantum numbers £j (cf. Ref. [0). 
In the tin nuclei, -Emax decreases from about 40MeV for the Si/2 states to zero for the 
kij/2 states. In Figs. ^l|-|33| different curves correspond to different cut-off energies, -Emax = 
-^'max + ^-E'max (^-E'max Varied bctwecn and 40MeV). Hence, Ai?max=40MeV corresponds 
to the cut-off energy -E^ax of 80MeV for the si/2 states, and 40MeV for the largest values 
of i. 

As expected, pairing energies depend significantly on the cut-off energy. Comparing 
results for Ai?max=40 MeV with those for Ai?max=0, one obtains differences of -Epair of the 
order of lO-i-20 MeV in the mid-shells. Due to the self-consistent readjustment of the p-h and 
p-p energies, the corresponding differences in the total energies (Fig. ^) are much smaller, 
2H-4MeV, but still significant. 

In the bottom panels of Figs. ^ and ^ are shown similar results for the renormalized 



strengths of the contact force (|2.2|). The values of Vq, quoted in the caption of Fig. have 
been obtained by requiring that the average neutron pairing gap in ^^°Sn, (Atv) =1.245 MeV, 
does not depend on AE'max- With such renormalized interactions, one obtains very small 
changes of total energies (Fig. bottom panel). The largest deviations do not exceed 
200 keV and 800 keV in stable and exotic isotopes, respectively, and can be safely disregarded 
when compared to all other uncertainties of methods used to extrapolate to unknown nuclei, 
or when studying the separation energies. 

Figure ^ shows the effective pairing-interaction strengths defined schematically as 
G'cfT=— (A)^/E'pair. The top panel presents the results obtained for Vo=— 160 MeVfm^ and 
for different values of Ai^max- One can see that the dependence on the cut-off energy is very 
weak, and the AE'max-dependence of (A)^ and -Epair cancels out in Gcs- (At A^=82 the pair- 
ing gap and the pairing energy both vanish, and hence the Geff values cannot be calculated.) 
A fixed value of Vq gives, therefore, a well-defined, cut-off-independent value of the effective 
pairing strength for every isotope. This result, together with the analysis of pairing gap 
distributions in Sec. |VI A| , demonstrates that calculations employing the volume contact p-p 
interaction are, in many respects, similar to those with the schematic seniority-pairing force 
(cf., however. Sec. [VF|) . 

The values of Gefr monotonically decrease with increasing neutron number N. The 



obtained dependence can be very well described by the simple Madland-Nix formula ||107 



G=llMeV/(ll+A^), while the Jensen-Miranda formula [|T08l, G=0.18 MeV[l-1.2/-2.8/2] 
{I={N—Z)/ A), gives a much faster decrease of Gcs with A^. (In a recent study ||109|] , based 
on a schematic finite-range force, the isospin dependence of G has been discussed. The 
authors found no sign of the P term suggested in Ref. ||108|] . This probably explains the 
disagreement seen in Fig. ^.) In both expressions we have normalized the multiplicative 
constants to obtain G=0.18MeV at N=Z =50. 
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The bottom panel of Fig. ^ shows similar results, but for the renormalized values of Vq, 
quoted in the caption of Fig. 0. The insert shows the values of Vq (dots) as function of 
Aii^max compared with the simple fit by the formula ( PT| ) with e„=ei+Ai?max, for e/=40 MeV, 
Co=430 MeVfm'^, and A=5.58MeV (solid line). One can see that the generic dependence 
of the renormalized values of Vq on Aii^max is fairly well reproduced, although the numerical 
constants Co and A obtained from the fit do not exactly correspond to the values inferred 
from the BCS theory with a constant density of single-particle states. 
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TABLES 



TABLE L Results of the HFB calculations with SkP force for the Si/2 neutrons in Sn. For 
the nth quasiparticle state, E^^^ and are the quasiparticle energy and the norm of the lower 
component, respectively. For the /ith canonical-basis state, v'^ is the occupation probability, and 



are, respectively, the average values of the p-h and p-p mean-field Hamiltonians [Eq. ( 4.22 )] 



and is the BCS-like quasiparticle energy defined in Eq. (4.24). All energies are in MeV. 
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FIGURES 



FIG. 1. Self-consistent spherical neutron densities PN{f) calculated with the SkP, SIII , and 
DIS interactions for selected tin isotopes across the /3-stability valley. 

FIG. 2. Self-consistent spherical neutron pairing densities PN{f) calculated with the SkP, 
SkP^, and DIS interactions for selected tin isotopes across the /9-stability valley. 

FIG. 3. Self-consistent spherical local neutron potentials UN{r) calculated with the SkP, SIII"^, 
and DIS interactions for selected tin isotopes across the /^-stability valley. 

FIG. 4. Self-consistent spherical local neutron pairing potentials UNir) calculated with the 
SkP and SkP*^ interactions for selected tin isotopes across the /3-stability valley. 

FIG. 5. The self-consistent HFB+SkP mass parameters M and M, and potentials U and U 
(central parts only), for neutrons in I20sn. 

FIG. 6. The HFB+SkP radial wave functions r(j)i{En,r) of the neutron 5^/2 single-quasiparticle 
states in i20Sn. Upper (i=l) and lower (i=2) components are plotted in the left and right columns, 
respectively. The numbers preceded by a times (x) sign indicate the scaling factors for the small 
wave function components. 

FIG. 7. The HFB+SkP radial canonical-basis wave functions r(/)^(r) of the neutron 5^/2 
single-particle states in i20sn. 

FIG. 8. Same as in Fig. I, but for the HF+BCS approach. 

FIG. 9. Weakly bound and unbound self-consistent single- neutron HF+SkP energies (^^j for 
^^'^Sn as function of -Rbox- Top and bottom panels show states of positive and negative parity, 
respectively. 

FIG. 10. Same as in Fig. ^, but for the canonical energies e„^j, Eq. ( f4.22a| ). 

FIG. 11. Self-consistent single-quasineutron HFB+SkP energies E^f^ (top panel) compared 
with the BCS-like canonical single-quasineutron energies E'^f[- (Eq. ( f4.24| )) (bottom panel) for the 
tin isotopes. 
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FIG. 12. The HFB+SkP spectral amplitudes ( ^ ) of the canonical states in i20sn 
with ^=1-5. The corresponding canonical energy is given in MeV and the occupation probability 
-u^ is displayed in parentheses. All the amplitudes for /li=1 have been assumed to be positive. 
This defines the relative phases of the spectral amplitudes for (shown by bars hashed in 

opposite directions). For E < —X the quasiparticle spectrum is discrete, while for E > —A it is 
represented by the discretized continuum. 

FIG. 13. Same as in Fig. but for ^^'^Sn. No discrete states appear for E < —X. 

FIG. 14. Same as in Fig. |l^, but for the //=l-3 canonical states. 



FIG. 15. Same as in Fig. y, but for ^^OSn 



FIG. 16. The self-consistent HFB+SkP single-neutron density PN{f) in ^^^Sn (in logarithmic 
scale) calculated with different values of i?box- The inset shows the same data in linear scale. The 
shaded line shows the asymptotic behavior given by Eq. ( ^.14a|) . 



FIG. 17. Same as in Fig. |T^, but for the neutron pairing density p]\[{r). The shaded line shows 
the asymptotic behavior given by Eq. (5.14b). 



FIG. 18. Self-consistent HFB-l-SkP neutron pairing potentials U{r) in ^^'^Sn (top panel) and 
^^^Sn (bottom panel) calculated with four values of -Rbox- The corresponding average gap values 
(A) [Eq. dJ)] are indicated. 

FIG. 19. Neutron rms radius r^v (top panel), average pairing gap (Ajv) (middle panel), and 
the Fermi energy Aat (bottom panel) calculated in the HFB-f-SkP model for ^^'^Sn (full circles) and 
^^^Sn (open circles) as functions of -Rbox- 

FIG. 20. Neutron densities Pn{t) (in logarithmic scale) calculated in the HFB+SkP and 
HFB+DIS models for ^^"^Sn, ^^^Sn, and ^^^Sn. 

FIG. 21. Same as Fig. ^ for the DIS interaction, but as a function of the number of oscillator 
shells Ns\^. For every A^'gh the oscillator-basis frequency cjq is adjusted as described in the text. 

FIG. 22. Same as in Fig. [l^, but for the HF+BCS+(A) approach with the constant pair- 
ing gaps as listed in Fig. |^ (top panel), the HF+BCS+SkP'' model (middle panel), and the 
HF+BCS+SkP model (bottom panel). In all these calculations the same pairing space (i.e., en- 
ergy cut-off) was used as in Fig. 
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FIG. 23. Average values of the neutron p-h and p-p potentials, enej and Antj [Eqs. ( [4.22|) ] in 
the canonical states calculated for ^^OSn in HFB+DIS (top), HFB+SkP'^ (middle), and HFB+SkP 
(bottom) model. Only the states with t;^^^>0.0001 are displayed. 



FIG. 24. Same as in Fig. |23|, but for ^^"Sn 



FIG. 25. Average neutron pairing gaps (Atv) [Eq. ( [4.6D ] calculated for SkP (solid), SIII'' 
(dashed), and DIS (dotted) interactions for the series of tin isotopes. 

FIG. 26. Sequences of nuclear single-particle levels for various potentials. Orbitals are labeled 
by the spherical quantum numbers. From left to right: (i) shell structure for a potential with 
spin-orbit term but with a very diffuse surface, (ii) the A'^osc=4 and 5 shells of the harmonic 
oscillator potential, (iii) no spin-orbit term, leading to a degenerate spin-orbit pattern as observed 
in, e.g., hypernuclei, and (iv) shell structure characteristic of nuclei near the stability valley. 

FIG. 27. Two-neutron separation energies for the A^=80, 82, 84, and 86 spherical even-even 
isotones calculated in the HFB-l-SkP model as a function of N/Z (lower scale, N=83) or Z (upper 
scale) . The arrows indicate the proximity of neutron and proton drip lines - see Fig. ^ for detailed 
predicted positions of two-particle drip lines. 

FIG. 28. Plot of the two-neutron separation energies S2n for all particle-bound even-even 
nuclei with A>16 and A^<208 calculated within the spherical HFB-|-SkP approach. Results for 
1905 nuclei are shown using the color codes spaced by 0.4 MeV for 52n<4MeV and by 2MeV for 
52n>4MeV. 

FIG. 29. Two- neutron separation energies S2n (top) and Fermi energies A at (bottom) for the 
Sn isotopes, calculated in the HFB approach with several Skyrme interactions and the Gogny-DlS 
interaction. 

FIG. 30. Pair transfer form factor, r^p(r), calculated directly from the HFB pairing density 
p{r) (solid lines), compared with the macroscopic form factor calculated from the derivative of the 
particle density Sp{r) (dashed lines). 

FIG. 31. Pairing energies -Epair in the tin isotopes calculated within the HFB-l-SkP'' model. 
Top panel shows the results for the fixed interaction strength Vo=— 160MeVfm'^ and for several 
cut-off energies AE^^ax added to the usual ^j-dependent cut-off energy -Emax [Hi- Bottom panel 
shows similar results when the values of Vq are renormalized to —158.64, —149.57, —145.41, and 
-142.01 MeV fm^ for AE,nax=W, 20, 30, and 40 MeV, respectively. 

FIG. 32. Same as in Fig. 0, but for the total energy relative to that obtained with A£'inax=0. 
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FIG. 33. Same as in Fig. 31, but for the effective pairing strength defined as Geff=— (A)^/£'pair. 
The insert shows the renormalized strength Vq compared with that given by Eq. ([Bl|). 
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